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ABSTRACT 

We combine a cosmological reionization simulation with box size of 100ft. _1 Mpc on a side and a 
Monte Carlo Lya radiative transfer code to model Lyman Alpha Emitters (LAEs) at z ~5.7. The 
model introduces Lya radiative transfer as the single factor for transforming the intrinsic Lya emission 
properties into the observed ones. Spatial diffusion of Lya photons from radiative transfer results 
in extended Lya emission and only the central part with high surface brightness can be observed. 
Because of radiative transfer, the appearance of LAEs depends on density and velocity structures in 
circumgalactic and intergalactic media as well as the viewing angle, which leads to a broad distribution 
of apparent (observed) Lya luminosity for a given intrinsic Lya luminosity. Radiative transfer also 
causes frequency diffusion of Lya photons. The resultant Lya line is asymmetric with a red tail. The 
peak of the Lya line shifts towards longer wavelength and the shift is anti-correlated with the apparent 
to intrinsic Lya luminosity ratio. The simple radiative transfer model provides a new framework for 
studying LAEs. It is able to explain an array of observed properties of z ^5.7 LAEs in Ouchi et al. 
(2008), producing Lya spectra, morphology, and apparent Lya luminosity function (LF) similar to 
those seen in observation. The broad distribution of apparent Lya luminosity at fixed UV luminosity 
provides a natural explanation for the observed UV LF, especially the turnover towards the low 
luminosity end. The model also reproduces the observed distribution of Lya equivalent width (EW) 
and explains the deficit of UV bright, high EW sources. Because of the broad distribution of the 
apparent to intrinsic Lya luminosity ratio, the model predicts effective duty cycles and Lya escape 
fractions for LAEs. 

Subject headings: cosmology: observations — galaxies: halos — galaxies: high-redshift — galax- 
ies: statistics — intergalactic medium — large-scale structure of universe — line: 
profiles — radiative transfer — scattering 



1. INTRODUCTION 



Mor e than four decades ago, iPartridge fc Peebles! 
(|1967ft proposed that prominent Lya emission re- 
processed from ionizing photons of young stars in 
galaxies can be used to detect high-redshift galax- 
ies. The first successful detections of high-redshift 
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imaging or with integral-field-units (IFU) spectroscopy. 
Owing to the high efficiency of target detection, LAEs 
naturally become objects for large surveys of high- 
redshift galaxies. Besides providing clues to the forma- 
tion and evolution of galaxies at the time when the uni- 
verse was still young, LAEs are an important tracer of 
the large-scale structure. The clustering of LAEs may be 
used to constrain cosmological parameters. In particular, 
the large- volume surveys such as the Hobby-E berly Tele- 
scope Dark Energy Experiment (HETDEX; IHill et afl 
2008) will enable the detection of the baryon acoustic 
oscillations (BAO) features (e.g... lEisenstein et al.ll2005l) 
in the LAE power spectrum. The BAO and the shape 
of the power spectrum can be used to measure the ex- 
pansion history of the universe at early epochs (z ~ 3), 
which constrains the evolution of dark energy and the 
curvature of the universe. 

LAEs are also a key probe of the high-redshift 
intergalactic medium (IGM), especially across the 
reionization epoch. The use of LAEs to learn about 
reionization has been the subj ect of intense study (e.g. , 
I Miralda-Escude fc Reei 119981: I Miralda-Escudl 119981: 
I Haiman fc Spaansl 119991: ISantos 20041: Haiman fc Cenl 
120051: IDiikstra et al.l 120071: IWvithe fc Cenl 12007ft . Suit- 
ably devised statistics, including luminosity function 
(LF) and correlation functions of LAEs, can be used to 
constrain the neutral fraction of the I G M during reion- 
izatio n (IMalhotra fc Rhoadsl 12004 IHaiman fc Cen 
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Mesinger fc Furlanettol 120081 llliev et all 120081: 
Daval et all I2008L 120091) . By comparing the LFs o f 
z ~ 5.7 and z ~ 6.5 LAEs, iMalhotra fc Rhoadsl (|200l 
conclude that r eionization was larg el y complete at 
z ~ 6 .5 (also see iDiikstra et all 120071 ). iMcQuinn et al.1 
(2007) show that, with the angular correlation function 
of the 58 available z ~ 6.6 L AEs in the Subaru Deep 
Field (|Kashikawa et all I2006D . limits may be placed 
on the IGM neutral fraction, favoring a fully ionized 
universe at z ~ 6.6. 

However, none of the previous work of LAEs men- 
tioned above used reionization simulations with concur- 
rent treatment of hydrodynamics plus radiative trans- 
fer of ionizing photons and Lya photons. Hydrody- 
namic and radiative transfer simulations provide realis- 
tic neutral gas distributions, and Lya radiative transfer 
yields detailed properties of the Lya emission. Realis- 
tic Lya radiative transfer calculations have been applied 
to high-redshift LAEs in cosmological simulations (e.g., 
iTasits iomi 200(|). The application, however, is limited to 
a few individual sources, which do not form a sample for 
st atistical study. 

IMcQuinn etafl (|2007t ) and llliev et all (|2008t ) studied a 
sample of LAEs in reionization simulations with cosmo- 
logical volume. However, the radiative transfer of Lya 
photons is treated in a simplistic way in their study: 
the observed Lya spectrum is modeled as the intrinsic 
line profile modified by exp(— r„), where r„ is the op- 
tical depth at frequency v along the line of sight. Al- 
though this exp(— Ty) model can yield insights into the 
properties of the observed Lya emission, such as the ef- 
fect of IGM on the observability of LAEs, it is far from 
a complete description of the radiative transfer of Lya 
photons. First, during the propagation, Lya photons ex- 
perience frequency diffusion, which is neglected by the 
simple exp(— t v ) model. The exp(— r„) model removes 
Lya photons at a given frequency according to the Lya 
optical depth, and no frequency change occurs for any 
Lya photon, therefore it does not yield correct Lya spec- 
tra. Second, the simple exp(— t v ) model does not account 
for the spatial diffusion of Lya photons either. LAEs 
in this model appear as point sources in Lya and there 
is no surface brightness information. Even if Lya pho- 
tons start from a point source, spatial diffusion due to 
radiative transfer would lead to an extended source. Ob- 
servationally, LAEs indeed appear to be extended and 
they are defined by a sur face brightness thr eshold in the 
narrow-band image (e.g.. lOuchi etalll2008h . Therefore, 
although the simple exp(— r„) model may provide use- 
ful insight, it likely falls short for predicting the detailed 
properties of the observed Lya emission from LAEs. 

To correctly understand high-redshift LAEs and use 
them for cosmological study, a full calculation of ra- 
diative transfer of Lya photons for a large sample of 
LAEs in cosmological reionization simulation is neces- 
sary, as will be evident later. In this work, we aim to 
perform detailed radiative transfer calc ulation of Lya 
photons (jZheng fc Miralda-Escudl 120021 ) fr om LAEs in 
a self-consistent fashion through radiation-hydrodynamic 
reionization simulations (jTrac et al.l 120081 ) . For this pa- 
per, we focus on studying statistical properties of z ~ 5.7 
LAEs and show how the radiative transfer calculation 
aids our understanding of the observed properties of 
LAEs. The clustering properties of LAEs from this study 



will be presented in another paper (Paper II; Zheng et 
al. in prep.). The paper is organized as follows. In § 2 
we review the cosmological reionization simulation used 
in this work and in § 3 we describe the Lya radiative 
transfer calculation. In § 4, we study in details the Lya 
emission from an individual source chosen from the sim- 
ulation box to gain a general view of the effect of Lya 
radiative transfer on the appearance of LAEs. Then, we 
present the statistical properties of LAEs in § 5, includ- 
ing their spectra and luminosity, from our modeling of 
an ensemble of sources in the simulation box. We com- 
pare our modeling results with observations for z ~ 5.7 
LAEs and discuss the implications in our understanding 
of LAEs. § 7 is devoted to identifying important physical 
factors in shaping the observed Lya emission of LAEs. 
We summarize and discuss the results in § 8. 

2. RADIATION HYDRODYNAMIC SIMULATION OF 
COSMOLOGICAL REIONIZATION 

In this work, we perform a Lya radiative transfer cal- 
culation to model LAEs. The sources and physical prop- 
erties of gas are taken from the outputs of a cosmological 
reionization simulation. 

The cosmological simulation (jTrac et al.ll2008l ) models 
cosmic reionization by using a hybrid approach to solve 
the couple d evolution of the dark matter, baryons, an d 
radiation (|Trac fc Pen! l200l [20061: ITrac fc Cenl [2007h . 
First, a high-resolution iV-body simulation was run, with 
3072 3 dark matter particles on a mesh of 11,520 3 cells in 
a box of 100/i~ 1 Mpc (comoving) on a side, and collapsed 
dark matter halos were identified on the fly. These ha- 
los are the sites to form sources of ionizing photons. The 
high resolution and large box size of the simulation make 
it possible to resolve small scale structures and to reduce 
sample variance for source statistics. 

Hydrodynamics and radiative transfer of ionizing pho- 
tons are simulated with moderate resolution (equal num- 
bers, 1536 3 , of dark matter particles, gas cells, and adap- 
tive rays). Within the limits of available computational 
resources, the multi-grid approach adopted in the simu- 
lation maximizes the resolution of the individual numer- 
ical components (gas and radiation) in order to model 
the corresponding physics adequately. The initial condi- 
tions are the same as in the high-resolution iV-body sim- 
ulation, and the high-resolution simulation is used only 
to generate a catalog of halos at each redshift step and 
to obtain the list of sources of ionizing radiation. These 
sources for the ionizing photons are then used in the lower 
resolution simulation. The sources are assumed to be 
Population II stars from starbursts ([Schaerei][2003) , and 
they are related to halos according to the prescription 
for sta r formation and emitted radiation in ITrac fc Cenl 
(2007). For each gas cell, the incident radiation flux is 
used to solve the temperature and ionization structure 
of each cell. For m ore details about the simulation, see 
ITrac fc Cenl (pOOl and ITrac et all (pOOSh . 

The simulation adopts a spatially flat ACDM cosmo- 
logical model with Gaussian initial density fluctuations, 
and the cosmological parameters are consistent with the 
Wilki nson Microwave Anisotropy Probe (WMAP) 5-year 
data (|Dunklev et all 12001 : O m = 0.28, Sl A = 0.72, 
n b = 0.046, h = 0.70, n s = 0.96, and a s = 0.82. 

3. Lya RADIATIVE TRANSFER CALCULATION 
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The outputs of the above simulation form the ba- 
sis for computing the radiative transfer of Lya pho- 
tons and studying LAEs. The radiative transfer of 
the resonan ce Lya line has been a subject of intense 
study (e.g... IHummerl [19621 I Auerl 119681: lAverv fe House! 
1968f IAdamslll972t iHarringtonll 19731 119741: lNeuieldl | 1990l 
199l[|Loeb fe Rvbickilll999|: lAhn et al.ll2000l 120011 200a 
Zheng fe Miralda-Escudel 120021: iDiikstra et al.l 120061 
Hansen fe Oh! l2006t iTasitsiomii [20061: iVerhamme et al.l 



20061: lLaursen et al I [20091 : iPierleoni et al.ll2009D . Owing 
to the complex nature of the geometry and gas distribu- 
tion in the cosmological realization we study, the Monte 
Carlo method of solving the Lya radiative transfer be- 
comes the na tural choice. We use the Monte Carlo code 
developed in iZheng fe Miralda-Escud c (2002), modified 
to use the simulation output, to solve the Lya radiative 
transfer in this study. This code has also been applied 
to study the fluorescent Ly a emission from the IG M in 
a hydrodynamic simulation (jKollmeier et al.ll2010l) . 
The code works as follows. 

(1) For each Lya photon, its initial position is drawn 
according to the emissivity distribution in the box 
(a superposition of point sources in the case pre- 
sented in this paper, see below). The initial fre- 
quency of the photon follows the Gaussian distribu- 
tion determined by the halo virial temperature (see 
below) in the rest frame of the fluid at the photon's 
position and its direction is randomly distributed. 

(2) An optical depth is then drawn from an exponential 
distribution. The spatial location along the chosen 
direction corresponding to this optical depth is de- 
termined from the distributions of neutral hydro- 
gen density, fluid velocity, and temperature along 
this direction. 

(3) At this location, the Lya photon encounters a scat- 
tering. The frequency and direction after the scat- 
tering are computed in the rest frame of the hydro- 
gen atom and then transferred back to the labora- 
tory frame. 

(4) With the new frequency and direction, (2)-(4) are 
repeated until the photon escapes from the system 
(see below). 

Lya photons are collected onto a three-dimensional (3D) 
array, which records the Lya spectra at each projected 
spatial location. At each scattering, as well as at the 
initialization, the possibility that the Lya photon escapes 
along the observational direction is computed and added 
into the array. In the end, the output array of the Monte 
Carlo Lya radiative transfer code forms an IFU-like data 
cube. Lya spectra (either ID or 2D) can be extracted 
from this data cube and Lya images can be obtained 
by collapsing the data cube along the spectral direction. 
For more details of the radiative transfer calculation , we 
refer the readers to IZheng fc Miralda-Escudel (|2002T ). 

In this paper, we focus on the simulation output at 
z = 5.7. The reionization is complete by that time in the 
simulation and the neutral hydrogen fraction is about 
9 x 10 -5 in the IGM. The neutral hydrogen density, tem- 
perature, and peculiar velocity fields in the simulation 
box are stored in a 768 3 grid, which feeds to the Lya 



radiative transfer code. The Hubble flow is added to 
the velocity field. LAEs are modeled to reside in dark 
matter halos. The positions and velocities of LAEs are 
from the halo list. To reduce source blending in the Lya 
image and spectra, Lya photons are collected with a spa- 
tial resolution finer than the above grid, with each grid 
resolved by 8 2 =64 pixels. The size of each pixel corre- 
sponds to 16.3/i _1 kpc (comoving) or 0.58". The resolu- 
tion of the 768 3 grid for gas properties used in the Lya 
radiative transfer calculation is a factor of two lower than 
in the hydrodynamical simulation, as a result of compu- 
tational efficiency consideration. The slight smoothing 
of gas fields may cause a smoothing effect in Lya surface 
brightness profile. However, since we use a finer grid to 
collect Lya photons and Lya sources are initially point 
sources (see below), the smoothing in Lya image is ex- 
pected to be much weaker than the smoothing in gas 
properties. 

The whole simulation box is divided into three layers 
along the line of sight, with the volume of each layer be- 
ing 100x100x33.33 /i" 3 Mpc 3 . The depth of each layer, 
33.33/i~ 1 Mpc, is close to the width of the narrow-band 
filter used to search for z ~5.7 LAEs in the 1 deg 2 field 
of the Subaru/ XMM-Newton Deep Survey (SXDS), and 
the area is almost identical to that of the survey as well 
(|Ouchi et al.H2008h . Therefore, we have three SXDS-like 
volumes at z *~5.7. For each layer, the output array on 
which the Lya photons are collected has a total range 
of 24A(rest frame) along the spectral direction, a width 
large enough to cover the Hubble expansion plus the pe- 
culiar velocities in the 33.33/i _1 Mpc width of each layer. 
The spectral resolution is set to be 0.lA(rest frame), cor- 
responding to 25km s . As a whole, the Lya radiative 
transfer results for each layer are saved in an array of 
dimension 6144x6144x240. 

We perform the Lya scattering calculation for all the 
halos above 5 x 10 9 h~ 1 MQ. Assuming ^2/3 of ionizing 
photons are converted to L ya pho tons (case-B re com- 
bination, lOsterbrockl |1989| ) and a iSalpeterl ()1955|) ini- 
tial mass function (IMF), the intrinsic Lya luminosity 
- ^intrinsic is related to t he star formation rate (SFR) as 
([Furlanetto et al.ll2005f ) 

(1) 



10 42 [SFR/(M o yr- 1 )] ergs" 



In the simulation, the result ant SFR under the adopted 
star formation prescription (jTrac fc Cenl 12007] ) is found 
to be tightly correlated with halo mass, 



SFR = 0.68[M h /(10 10 /i" 1 Af o )] M Q yr- 



(2) 



So the intrinsic Lya luminosity and halo mass are al- 
most interchangeable in our model and in our descrip- 
tions of the results. The relation in equation ([2]) holds 
at z ^5.7, and there is a redshift dependence (see 
iTrac fc Cenl [20071) . Since the ultra-violet (UV) luminos- 
ity is also proportional to SFR, the halo mass to UV (or 
intrinsic Lya ) light ratio is approximately constant in 
our model. 

For each halo, Lya photons are launched at the halo 
center. The point source assumption is reasonable. Lya 
emission origi nates from reprocessed ioni zing photons of 
massive stars ([Partridge fc Peebles! [1961 . The ionizing 
photons ionize the neutral hydrogen atoms in the in- 
terstellar medium (ISM) and the case-B recombination 
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Fig. 1. — Lya image and spectra of a single z ~ 5.7 LAE, randomly chosen in the simulation box. Left: The Lya surface brightness 
distribution after the radiative transfer. The physical size of the region is indicated by the scale bar and the separation of the two adjacent 
ticks on the axes correspond to 50kpc. Two circles delineate the apertures for the spectra in the right panels. The virial radius of the host 
halo (~ 10 11 /! - 1 Mq) is about 26kpc, slightly larger than that of the inner circle. Right: Lya spectra (solid curves) at different radii to the 
source center. Panel (a) shows the spectrum within the smaller circular aperture shown in the image, and Panel (ft) for photons in between 
the smaller and the larger circles and Panel (c) for photons outside of the larger circle. The dotted curve in each panel shows the intrinsic 
Lya line profile, which would be observed if there were no scatterings. The spectra are arbitrarily normalized, since we concentrate on the 
profiles. The wavelength is shown as the difference to (1 + z)Aq, where z = 5.7 and Ao=1216Ais the rest-frame wavelength of Lya . 



has a probability of ~2/3 of ending up as Lya photons 
(Ostcrbrock 1989). We aim to solve the radiative trans- 
fer in the circumgalactic and intergalactic environments, 
and the initial Lya photons launched in our model cor- 
respond to photons just escaping from the ISM whose 
spatial distribution closely follows the UV light of galax- 
ies. From HST/ACS observations of z = 5.7 LAEs, 
iTaniguchi et al.l (j2009T ) find that in the broad-band (rest- 
frame UV) images, LAEs are compact sources with sizes 
of less than 1 kpc, smaller than the pixel size in our mod- 
eling. Therefore, our assumption of a point source for the 
initial Lya emission is justified. 

The initial frequency of the Lya photons in the rest 
frame of the halos is assumed to follow a Gaussian 
distribution with the width corresponding to the virial 
temperature of th e halo, T V1I = GAf/ l /im///(3fci? V ir) 
(|Trac fc Cenl 12007ft . where [i ~ 0.59 is the mean 
molecular weight. The rms line width is CTi n ;t = 
31.9[M ?i /(lO lo /i- 1 M )] 1 / 3 kms- 1 . This width is about 
half of the circular velocity at the virial radius. This 
is a rather conservative assumption and we will test 
and discuss the effect of increasing the initial width on 
our results. The total number of Lya photons drawn 
for each halo is N 1 =max{SFR/(M Q yr- 1 ),10 3 }, and 
these photons are given a weighting factor w to convert 
the photon number to the Lya luminosity of the halo, 

-^intrinsic — wN*y . 

To fully account for the effect of the IGM on Lya radia- 
tive transfer, we impose the periodic boundary conditions 
of the simulation in our Lya radiative transfer calcula- 
tion. For each Lya photon, we stop the scattering calcu- 
lation when it reaches a distance of half of the box size 



(L = 100/i _1 Mpc) from the initial position in any of the 
three principle directions. At this distance, the Hubble 
expansion leads to a fractional shift in Lya wavelength 
of the order of AA/A = 0.5H(z)L/(l + z)/c ~ 1.6 x 10~ 2 , 
corresponding to a velocity of ^5000kms _1 . This is 
much larger than typical values of peculiar velocity of 
halos and the shift caused by frequency diffusion, and is 
more than sufficient to ensure that the photon will no 
longer interact at the Lya line. In fact, most of the time 
the photon is last scattered at a distance a few comoving 
Mpc away from the source (e.g., see Fig. [T]). 

4. DETAILED STUDY OF AN INDIVIDUAL LAE 

Before presenting Lya radiative transfer and statistical 
results for all the LAE sources in the whole simulation 
box, we first examine in detail the radiative transfer re- 
sults for an individual source to aid our understanding 
of the general effects of Lya scattering. 

We randomly chose a halo in the simulation box, which 
has a mass of ~ 10 u ft. _1 M Q . The size of the virialized 
halo is about 26kpc, slightly larger than the inner cir- 
cle in the left panel of Figure [1] As mentioned in § |3l 
Lya photons are assumed to initially start from a point 
source, located at the halo center. The left panel of Fig- 
ure [1] shows the Lya surface brightness distribution after 
Lya scattering from this source. Because of the radiative 
transfer, the initial point source becomes an extended 
source and a roughly spherical halo of scattered Lya pho- 
tons emerges. This scattered Lya halo is similar to the 
o ne around a poin t sourc e befor e reionization described 
in lLoeb fc Rvbickl (fl999l) . While ILoeb fc Rvbickil ((19991) 
assume a uniform, zero temperature IGM undergoing 
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Fig. 2. — Lya images (left) and spectra (right) of a single z ~ 5.7 LAE observed along two opposite directions. Mirroring reflection has 
been applied so that the two images have the same orientation, fn the image panels, the physical size of the region is indicated by the scale 
bar and the separation of the two adjacent ticks on the axes correspond to 50kpc. In each of the right panel, the dotted curve shows the 
intrinsic Lya line profile, which would be observed if there were no scatterings. The solid black curve shows the spectrum of all scattered 
photons and the solid red curve is the spectrum near the center (within the circular aperture shown in the image panel). All quantities are 
in the observer's frame. The wavelength is shown as the difference to (1 + z)Ao, where z = 5.7 and Ao=1216Ais the rest-frame wavelength 
of Lya . The relative shift of wavelength ranges in the top and bottom spectra panels reflects the change in the viewing direction, which 
leads to the change in the source position with respect to the box center and that in the velocity direction with respect to the observer. 



Hubble expansion, we use a realistic distribution of gas 
density, temperature, and velocity around a star-forming 
halo. This causes deviations from spherical symmetry in 
the surface brightness profile. The scattered Lya surface 
brightness drops as the radius increases. The blue-green 
extended Lya emission seen in Figure[T]is out in the IGM, 
corresponding to the scattered photons as they travel to 
the region where the Hubble expansion compensates for 
their blueshift acquired by the scatterings in the infall 
region of the halo. The sharp edge around a radius of 
~0.5Mpc reflects the frequency of the "bluest" photons 
coming out of the central (infall) region before encounter- 
ing the IGM. The bluer the photons are, the farther they 
can travel in the IGM before redshifted to the line center 
and significantly scattered. In practice, we can only ob- 



serve the very inner part of the extended Lya radiation, 
where the surface brightness is high. The extended Lya 
halo would merge with those from neighboring sources, 
forming a Lya background. We will describe how we 
identify sources in the simulation in § 15.11 

In the right panels of Figure [TJ we show Lya spec- 
tra at different radii of the source. The dotted curve in 
each panel is the intrinsic Lya line profile, assumed to be 
Gaussian with the width determined by the virial tem- 
perature of the dark matter halo hosting the source. It 
would be the observed spectrum if Lya photons streamed 
out of the source without any scatterings. Note that 
the wavelength shown in the plot is the difference from 
(1 + z)Xq, where z = 5.7 and Ao=1216Ais the rest-frame 
wavelength of Lya . The offset of the peak of the initial 
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Fig. 3. — Same as Fig.[2]but with the peculiar velocity field turned off. 



line profile from zero is caused by a combination of the 
Hubble velocity with respect to the center of the simu- 
lation box and the peculiar velocity of the source. On 
average, Lya photons are first scattered by neutral hy- 
drogen atoms in the infall region around the LAE host 
halo (see § [7] and Fig. [20]) . The inner infall region has 
an inverted Hubble-like contraction velocity distribution 
and Lya photons escaping at the radius of maximum 
infall velocity have their frequ ency most likely shifted to 
the b lue side of the line center (|Zheng fc M iralda- Escudel 
2002). Then Lya photons experience scatterings in the 
region with decreasing infall velocity (outer infall region) 
and finally hit the Hubble flow. The subsequent scat- 
terings on average shift the frequency of Lya photons 
redward. 

Most of the observed Lya photons from regions with 
high surface brightness near the center of the source shift 
to the red side (Fig. [l}z) with respect to the intrinsic dis- 
tribution. These photons are most likely to have had for- 
ward scatterings along the line of sight. The outward in- 



creasing gas velocity (from the outer infall region and the 
Hubble flow) makes it easier for photons that have red- 
ward frequency shifts to escape. At larger lateral radii, 
(Fig. QJ and Fig. [T];), we have more contributions from 
photons that travel perpendicular to the line of sight but 
are scattered to the line-of-sight direction. The frequency 
after scattering would be near the incoming frequency 
seen in the rest-frame of the atom at the scattering ra- 
dius. Since most of the scatterings would happen around 
the radius where the Hubble expansion velocity redshifts 
the Lya photons to the line-center frequency and the 
line-of-sight neutral column density become smaller at 
larger radii, the observed photons at large projected radii 
would appear bluer than those observed near the center, 
as seen in the spectra at large radii. 

Only the very central part of the extended Lya ra- 
diation can be practically observed as the LAE source. 
The transmitted Lya flux depends on the gas distribu- 
tion and kinematics in the halo vicinity, and the view- 
ing angle. We perform a few tests to see the influence 
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of the gas properties in the observed flux of the central 
part. Figure [2] compares Lya images and spectra of the 
above source observed in two opposite directions. Mir- 
roring reflection has been applied to one image so that 
the two images have the same orientation. While the two 
images have similar spatial extent, there is a large differ- 
ence (about a factor of 7) in the flux inside the central 
aperture. The difference can be clearly seen from the 
spectra. The red curve in each of the right panels is the 
spectrum extracted from the circular aperture near the 
source center. 

Evidently, the differences in the intervening gas distri- 
bution, namely the neutral hydrogen density and pecu- 
liar velocity distributions, have a dramatic effect on the 
observed Lya flux. To test this we perform a scattering 
calculation with the peculiar velocities of the source and 
neutral gas set to zero, while keeping the Hubble expan- 
sion. Figure [3] compares the resultant images and spectra 
observed in the two opposite directions. Any remaining 
differences between the results of different viewing direc- 
tions should be caused only by the anisotropic density 
or temperature distribution around the source. We see 
that, with the peculiar velocity field turned off, the dif- 
ference in the fluxes from the central aperture between 
the two lines of sight becomes much smaller, a factor of 
~ 1.5 (versus ~ 7 in the case with peculiar velocity). The 
spectra from the central aperture also look more similar 
to each other. 

With the peculiar velocity turned off, the surface 
brightness profile of scattered Lya photons appears to 
be more concentrated than that in Figure [3] This is 
mainly a consequence of the disappearance of the infall 
region around the source by artificially setting the pe- 
culiar velocity to zero. If the peculiar velocity is not 
switched off, Lya photons climbing out of the infall re- 
gion (before reaching the IGM dominated by the Hub- 
ble expansion) would on averag e have shifted blueward 
(|Zheng fc Miralda-Escudell200a ) with respect to the line 
center. Compared to the case without the blueward shift 
(e.g., the initial Gaussian profile when the peculiar ve- 
locity is turned off), these bluer photons would travel a 
larger distance in the IGM before redshifting back to the 
line center and experiencing strong scatterings. There- 
fore, we see a more extended Lya emission in the case 
with the peculiar velocity. 

The above tests show that peculiar velocity plays an 
important role in the scattered Lya brightness profile 
and the transmitted flux near the source center. In § [?1 
in additional to the peculiar velocity, we identify other 
factors in affecting the Lya transmission and statistically 
study their role in the observability of LAEs. 

The results from the individual LAE source show that, 
as a consequence of radiative transfer, Lya photons ex- 
perience both spatial diffusion and frequency diffusion. 
An intrinsic point source of Lya emission appears ex- 
tended and the Lya spectra differ substantially from the 
intrinsic Gaussian profile. The spectra from the central 
aperture, which is the part that is most observable, do 
not have a simple and clear relation to the initial line 
profile, owing to the frequency sh ift caused by scatter- 
ings. In some p revious work (e.g.. iMcQuinn et alj|2007t 
llliev et aT1l2008D , the observed Lya spectrum is modeled 
as the intrinsic Gaussian profile modified by exp(— t„) 
with t„ being the optical depth at frequency v along 



the line of sight. Such a simple model does not account 
for the frequency and spatial diffusion of Lya photons 
caused by radiative transfer. The resultant line profile in 
this simple model looks like a truncated Gaussian pro- 
file, with only the red tail transmitted. Our detailed Lya 
radiative transfer, on the other hand, shows that the ob- 
served Lya line profile near the source center is more 
complicated and the redward frequency shift is more than 
that in the simple treatment. The simple radiative trans- 
fer model may yield trends in some results that are quali- 
tatively i n accord wi t h deta iled transfer calculations. For 
example, llliev et all ()2008l ) also find that peculiar veloc- 
ity is important in determining the observability of LAEs. 
However, as we show in § El the lack of frequency and 
spatial diffusion in the simple model means that it cannot 
capture the full picture of Lya emission from LAEs for 
detailed prediction and understanding of observed Lya 
features. 

5. STATISTICAL PROPERTIES OF Lya SPECTRA AND 
LUMINOSITY OF LAES 

We perform Lya scattering calculation for all the 
sources residing in halos above 5 x 10 9 /i _1 M Q in the 
whole (lOO/i Mpc) box. In this section, we describe 
how we identify LAEs from the post-scattering outputs 
and study their statistical properties. 

5.1. Source Identification 

Figure |4] shows the Lya image for sources in 
one-third of the simulation box. It has an 
area of 100/i _1 Mpcxl00/i _1 Mpc and a thickness of 
33.33/i Mpc. The slice matches the sky coverage (1 
deg 2 ) of the SXDS and the depth corresponds to the 
width of the narrow-band filter (AA=120A; lOuchi et al.l 
2008) for z = 5.7 LAEs. Therefore, the image can be 
regarded as an idealized, continuum-subtracted narrow 
band image of the z = 5.7 LAEs for SXDS-like sky cov- 
erage and depth. From the whole simulation box, we 
have three realizations of such a survey. The periodic 
boundary condition of the simulation is imposed in our 
modeling, which can be clearly seen in Figure [U 

Because of Lya radiative transfer, LAEs are no longer 
point sources in our model. We need to find a way 
to define the sources in order to study their statistical 
properties. Our identification of sources is motivated by 
the procedure used in detecting LAEs in real observa- 
tions. For z = 5.7 LAEs in the SXDS, a threshold sur- 
face brightness of 2.64 x 10 _18 ergs _1 cm~ 2 arcsec -2 in 
the narrow-band image (including continuum) is adopted 
for detecting them (M. Ouchi, private communication). 
LAEs are identified by grouping pixels above this thresh- 
old. The observed rest-frame Lya equivalent width dis- 
tribution of z ~ 5.7 LAEs peaks around 60Aand i s 
skewed to large values (Fig. 23 in lOuchi et al.1 [2008). 
Since the rest-frame width of the narrow-band filter is 
120A/(1 + z), the continuum contribution to the sur- 
face brightness are likely to be less than 30%. Our 
model does not include the continuum component. In 
principle, we could model the continuum based on the 
star formation history in the simulation, but the cor- 
rection is small and it is not the main uncertainty of 
our model (as shown later in this paper). Therefore, 
we simply make a correction of 1/3 to remove the con- 
tinuum contribution to the threshold surface bright- 
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Fig. 4. — Scattered Lya emission from LAEs. The Lya surface brightness distribution is shown for sources in one third of the simulation 
box, with 100ft - ^^Mpc on a side and a thickness o f 33.33fe~ 1 Mpc. T he area matches that of the SXDS and the depth corresponds to the 
width of the narrow-band filter for z ~ 5.7 LAEs (Ouchi ct al. 2008). The morphology and spectra of three sources in the zoom-in region 
are shown. 



ness and adopt 1.80 x 10 18 ergs 
continuum-subtracted surface brightness for detecting 
SXDS LAEs. 

As shown in § El there is significant uncertainty 
in modeling Lya luminosity. To be conservative, 
we set a lower threshold surface brightness, 1.8 x 
10 _19 ergs~ 1 cm~ 2 arcsec~ 2 , for identifying LAEs in our 
model. To have a full picture of the observability of LAEs 
(e.g., comparison between intrinsic and observed Lya 
emission), we also include the halo position information 
in the source identification. For a source with the pro- 
jected position known from the halo catalog, we first find 
the corresponding pixel in the Lya image. Starting from 
this root pixel, we then link the surrounding pixels with 
surface brightness above the threshold through a Friends- 
of-Friends (FoF) algorithm with the linking length equal 
to the size of the pixel. The directions to link the pix- 
els are only horizontal and vertical in the image. An 
LAE source is then defined by all the linked pixels, and 
the spectra at each pixel's position are added together to 
form the spectra of the source. In the case that the root 
pixel and its surrounding pixels all have surface bright- 



nesses lower than the threshold, the flux and spectra from 
this root pixel are adopted for the source. As with other 
applications of the FoF algorithm, there are chances that 
two individual sources are bridged together. In this case, 
we also assign the flux and spectra of the root pixel to 
the intrinsically fainter source. Such cases are rare and 
the correction does not affect any of our statistical study. 

Figure[5]shows Lya images and spectra of a few z ~ 5.7 
LAEs from our model. The host halo masses of these 
sources are above 3 x 10 10 h~ 1 M ( ? ) . The corresponding in- 
trinsic Lya luminosities are above 2 x 10 42 ergs _1 (equa- 
tions [Tj and [5] ) , roughly in the luminosity range probed 
by current LAE surveys like SXDS. From top to bottom 
panels and left to right panels, they are arranged in or- 
der of decreasing intrinsic Lya luminosity (halo mass). 
Most of the sources appear to be roughly round with faint 
substructures around them, which are a combination of 
reflected Lya emission by clumps/filaments of neutral 
gas or Lya emission from fainter sources. The sizes and 
morphologies of the LAEs in our model are remarkably 
similar to those in the narro w-band images of z ~ 5.7 
LAEs in SXDS (e.g., Fig.5 of IQuchi et alJI2008t narrow- 
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Fig. 5. — Images and spectra of a few z ~ 5.7 LAE sources in our model. In each of the panels of spectra, the dotted curve is the intrinsic 
line profile, which would be observed if Lya photons escaped without scattering. The black solid curve is the spectrum in our model with 
the full Lyc* radiative transfer. For comparison, the solid red curve is the spectrum with a simple treatment of radiative transfer, which 
modifies the intrinsic profile by multiplying exp(— t„) with r v the line-of-sight optical depth. For an easy comparison across panels, in each 
panel, all three spectra have the same constant horizontal shift so that the line center of the intrinsic profile is at zero. In each spectra 
panel, the label starting with "M" denotes \og(M h ), where the host halo mass M h is in units of h~ Mq. See text for details. 



band images in Fig.3 of lTaniguchi"eTaLl [20091 ) . 

Solid curves in the spectrum panels in Figure [5] are 
the corresponding spectra for the shown LAEs. In each 
spectrum, the Lya line is clearly asymmetric, skewed to- 
wards the red. The line profiles resemble the observed 
ones for the SXDS z ~ 5.7 LAEs (Fig.5 of lOuchi et all 
2008). However, the observed Lya lines appear to be 
much broader with less sharp blue edges. The differ- 
ence can be simply attributed to the spectral resolution: 
in the observer's frame, the observation typically has a 
resolution 8-15A (|Ouchi et aT1l2008l ). while the resolution 
for our modeled spectra is 0.67A. 

The Lya lines with the full radiative transfer show a 
clear distinction from the lines with a simple treatment of 
the radiative transfer, namely the exp(— t v ) model. For 
each source, the red curve is the Lya spectrum from the 



exp(— r„) model, which is essentially the intrinsic Gaus- 
sian profile truncated below a certain wavelength. Al- 
though it displays a similar asymmetry, the flux is usu- 
ally significantly lower than that with the full radiative 
transfer. Importantly, the Lya line from the full radiative 
transfer model has a larger redward shift than that in the 
simple model, an effect that can only be properly mod- 
eled with detailed radiative transfer calculations. This 
is primarily because detailed radiative transfer leads to 
frequency diffusion, causing some of the original photons 
closer to the line center to diffuse out to the wings. The 
shift in frequency not only results in smaller scattering 
optical depth but also is accompanied by spatial diffu- 
sion, both leading to larger transmitted flux near the 
center. 
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Fig. 6. — The apparent shift in the LAE Lya line peak with respect to the intrinsic one. Because of radiative transfer, the wavelength 
^pcak a t the peak flux of the observed Lya line is not A c that corresponds to the center of the intrinsic (Gaussian) line profile. The 
apparent shift is defined as the wavelength difference. The top axis of the left panel marks the shift in comoving unit, which is of the order 
of fc _1 Mpc. Left: The distribution of the apparent shift of the peak wavelength as a function of LAE host halo mass (in units of h~ 1 Mq). 
The apparent wavelength shift can translate to an apparent shift in the redshift-space source position AZ along the line of sight, which is 
labeled in the top axis. Right: Similar to the left panel, but the peak wavelength shift is in units of the intrinsic line width a\, which is set 
by the virial temperature of the host halo in our model. 
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Fig. 7. — Similar to the left panel of Fig. [6j but for the appar- 
ent shift in the LAE Lya line peak as a function of the observed 
(apparent) Lya luminosity (in units of ergs -1 ). 

5.2. Shift in the Peak of Lya Spectra 

Together with the assumed intrinsic properties, LAEs 
identified in the post-scattering IFU-like data cube from 
our model enable a statistical study of the relations be- 
tween the observed and intrinsic quantities, which in- 
clude spectral features and luminosity. 

As shown in Figure^ owing to radiative transfer effect, 
the peak in the observed Lya spectra is at a wavelength 
longer than that in the intrinsic spectra. In the left panel 
of Figure^ we plot the distribution of the shift as a func- 
tion of host halo mass. The broad distribution of the 
shift reflects the distribution of circumgalactic and inter- 
galactic environments (density and velocity structures), 



which affect the radiative transfer of Lya photons (see 
§[7J). This apparent wavelength shift in the peak depends 
on host halo mass, or intrinsic luminosity of the source 
given that star formation rate is tightly correlated with 
halo mass in the reionization simulation. In general, the 
distribution is skewed to large shifts. For sources in lower 
mass halos, the distribution is narrower and the average 
shift is smaller. For sources in halos of ~ 10 10 }i~ 1 Mq, 
the median shift corresponds to a velocity of ~ 70km s _1 , 
while for - 1O 11 /i~ 1 M halos, the value is - 200kms~ 1 . 
If the observed Lya emission was from photo ns backscat- 
tered from the far side of gal actic wind (e.g,. IFranx et al.1 
U997tlAdelberger et al.H2003f ). the above shift would lead 
to an overestimate of the wind velocity as long as it is es- 
timated by the apparent velocity difference between the 
Lya emission and optical emission/absorption lines. The 
apparent wavelength shift also translates to a shift in the 
apparent position/redshift of the source along the line of 
sight (top axis of the left panel of Fig. [5]). This position 
shift and its scatter would result in slight distortion and 
smoothing in the clustering of LAEs in redshift space 
(see Paper II for more details). 

In the right panel of Figure [6l the apparent shift of 
the Lya line peak is put in units of the intrinsic Lya line 
width <7\ , which is assumed to be determined by the halo 
virial temperature. Although the trends seen in the left 
panel are still evident, the variation of the distribution 
as a function of halo mass becomes weaker. Roughly 
speaking, the median shift is about 3o~\ and the scatter 
is about 0.6-0.9<7a- In terms of velocity, the median shift 
is approximately 100[M h /(10 10 ft.~ 1 M Q )] 1 / 3 kms" 1 . 

The distribution of Lya peak shift we discuss so far is 
as a function of intrinsic Lya luminosity (i.e., halo mass). 
From the point of view of observation, it is of great in- 
terest to show the distribution as a function of the ob- 
served Lya luminosity. Figure [7] plots the distribution of 
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Fig. 8. — Relation between the observed (apparent) and the intrinsic Lya luminosities of LAEs. Panel (a): The joint distribution of 
apparent and intrinsic Lya luminosities. Adjacent contours differ by a factor of 2 in contour levels. A vertical cut in this plot gives the 
probability distribution function (PDF) of apparent luminosity at a given intrinsic luminosity [Panel (b)], and a horizontal cut gives the 
distribution of intrinsic luminosity at a given apparent luminosity [Panel (c)]. Only halos above 5x 10 9 h~ 1 MQ are considered in our model, 
which corresponds to intrinsic luminosity above ~ 10 41,5 ergs _1 . As a results, the apparent luminosity is complete above ~ 10 41,2 erg s — 1 . 



Lya peak shift as a function of observed (apparent) Lya 
luminosity. As shown later in § 15.31 at fixed intrinsic 
Lya luminosity, the observed (apparent) Lya luminos- 
ity has a broad distribution, and vice versa. The peak 
shift distribution at fixed observed luminosity is there- 
fore contributed by sources residing in halos of a broad 
range of mass. We note that the model luminosity in the 
plot should be increased by about 0.7 dex to match the 
z ~ 5.7 observation (see §|S]). 

5.3. Apparent and Intrinsic Lya Luminosities 

The observed Lya flux F^ ya from an LAE source in the 
narrow-band image comes from the central part, where 
the surface brightness is high. That is, only a fraction of 
the extended radiation composed of scattered Lya pho- 
tons can be observed. The observationally inferred Lya 
luminosity ^apparent = 47T_D| / FL ya is therefore expected 



to be lower than the intrinsic Lya luminosity ^intrinsic, 
where Dl is the luminosity distance and isotropic emis- 
sion is assumed in computing the apparent luminosity 
from the observed flux. 

We compare the intrinsic and apparent Lya luminosi- 
ties in our model. Figure EJi gives the joint distribution 
of I/intrinsic and Iz a pparent- From the joint distribution, we 
can obtain the distribution of ^apparent at a fixed ^intrinsic 
or vice versa, through a vertical or horizontal cut. Since 
we only model LAEs in halos above 5 x 10 9 /i _1 Mq , we are 
limited to sources with ^intrinsic above ~ 10 41 ^ergs _1 . 
When considering the observed luminosity, sources are 
complete for L a pparcnt > 10 41 - 2 ergs _1 . We note that the 
Lya luminosity limits may change if the assumed IMF 
and SFR differ from those in our model. 

In the luminosity range probed by our model, the ap- 
parent Lya luminosity peaks at a few percent of the in- 
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Fig. 9. — Correlation between apparent to intrinsic luminosity 
ratio and wavelength shift at peak Lyct flux. The peak wavelength 
shift is in units of the intrinsic line width cr\, which is set by the 
virial temperature of the host halo in our model. Contours show 
the correlation and they enclose the 68% and 95% distribution, 
respectively. The solid curve is the median wavelength shift as a 
function of luminosity ratio. The two dotted curves delineate the 
upper and the lower quartilcs. 

trinsic one and is broadly distributed (Fig(5f0 . The scat- 
ter reflects differences in neutral gas distributions (den- 
sity, velocity, and temperature) around sources of the 
same intrinsic luminosity (§0. The distribution shifts 
slightly towards higher values for sources of lower intrin- 
sic luminosity. If we define the ratio of the apparent to 
intrinsic luminosity as flux suppression, the suppression 
on average appears to be smaller for intrinsically fainter 
sources, which seems counter-intuitive. Such a shift is 
a consequence that the environment of low mass halos 
are on average less dense than that of massive halos, and 
that the environment is important in shaping the observ- 
ability, to be discussed in detail in § [7] 

From the point of view of observation, it is interesting 
to ask what the observed Lya luminosity implies about 
the intrinsic one. Figure \EJp shows the intrinsic lumi- 
nosity distribution at a given apparent luminosity. The 
distributions for different values of £ apparen t are similar 
in terms of the intrinsic to apparent luminosity ratio. In 
general, the intrinsic luminosity is about 3-12 times the 
observed luminosity. In other words, a large fraction of 
the escaped Lya photons are invisible. For estimating 
Lya escape fraction from observations, this is a system- 
atic factor that needs to be taken into account. 

We find that the flux suppression is correlated with 
the shift in the peak of Lya profile (§ 15. 2[) , as shown in 
FigO For sources with a larger suppression in Lya flux, 
the peak of the spectra shifts more towards red. This 
correlation has only a weak dependence on halo mass or 
intrinsic luminosity and in Figure [9] all sources in our 
model are included. Clearly, the correlation is a con- 
sequence of the radiative transfer: Lya photons diffuse 
more in frequency as they experience more scatterings. 
The correlation is driven by the dependence of the Lya 
radiative transfer on environments, i.e., the circumgalac- 
tic and intergalactic density and velocity structures (see 
§0- 



It is worth pointing out that the simple exp(— t v ) 
model can only give qualitatively trends seen in our re- 
sults. In Figure [lOl we compare the L appai cnt-£intrinsic 
distribution from our model of full radiative transfer with 
those from the exp(— r„) model. The exp(— r„) model in 
Figure fTOh assumes the same intrinsic Lya line width as 
in our model, which is set by halo virial temperature. 
It is evident that at the same intrinsic Lya luminos- 
ity, the exp ( — t v ) model leads to much lower apparent 
Lya luminosities than the full calculation. In particular, 
the suppression from the exp(— r„) model becomes much 
stronger for sources of higher intrinsic Lya luminosity 
(or halo mass) because of the high de nsity and peculiar 
velocity. The trend is similar to what llliev et all (2008) 
find. As they point out, peculiar velocity plays an im- 
portant role in shaping the luminous end of the observed 
Lya luminosity function. The frequency and spatial dif- 
fusions in the full calculation can compensate the density 
and peculiar velocity effect, weakening the suppression. 
Compared to the exp(— t„) model, the suppression from 
the full calculation does not become much stronger for 
sources of higher intrinsic Lya luminosity (also see Fig- 
ure HU). 

The exp(— r„) model in Figure [TOb adopts an intrinsic 
Lya line width 2.3 times that used in Figure [TUb . which 
corresponds to the circular velocity at halo virial radius. 
This value of int rinsic line width is used in some pre- 
vious work (e.g., iMcQuinn et all 120071 : llliev et all [20081 ) 
with the exp ( — t v ) model. This model boosts the appar- 
ent Lya luminosity, but the suppression is still stronger 
than the full calculation at high halo mass end. In addi- 
tion, the distribution of apparent Lya luminosity at fixed 
intrinsic Lya luminosity is much narrower than that from 
the full calculation. We do not have results from a full ra- 
diative transfer calculation with the larger intrinsic Lya 
line width yet, but we expect that the difference between 
such a full calculation and the exp(— r„) model is similar 
to that seen in Figure [TUb . We caution that Figure [TOb 
does not show an apple-to-apple comparison, since the 
exp(— r„) model and the full calculation assume different 
intrinsic line widths. Nevertheless, it indicates that mod- 
ifying the exp ( — t u ) by varying the intrinsic line width 
does not lead to a match to the full radiative transfer 
calculation. 

To further see the difference between the full radiative 
transfer model and the exp(— t„) model with the same 
intrinsic Lya line width setup, we compare the apparent- 
to-intrinsic Lya luminosity ratios predicted from the two 
models on a one-to-one basis (Figure [TTJ) . We summarize 
the scatter plot by showing the median ratio (solid curve) 
from the exp(— t v ) model as a function of the ratio from 
our model (the "true" ratio) , together with the lower and 
upper quartiles (dotted curves). Thin and thick curves 
are for sources in halos below and above 10 10 5 h~ 1 Mo, 
respectively. In general, there is a trend that the ra- 
tio from the exp(— r„) model increases with the "true" 
value and this trend seems to break down in massive ha- 
los at high values of the "true" ratio. In a limited range 
(around L appar cnt/£intrinsic ~ 10~ 2 ), the median ratio 
from the exp(— t„) appears to be a constant shift from 
the "true" ratio. However, even if we apply a correc- 
tion to account for the shift, the exp(— r„) model would 
underpredict the luminosity ratio outside of the above 
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Fig. 10. — Comparison of the relations between the observed (apparent) and the intrinsic Lya luminosities from different models. The 
dotted contours are from our model with full calculation of Lya radiative transfer (the same contours as in Fig. [St). The solid contours in 
panel (a) are from the simple exp(— t v ) model with the initial Lya line width given by halo virial temperature, as adopted in our model 
with full radiative transfer calculation. The solid contours in panel (b) are also from the exp(— t„) model, but the initial Lya line width is 
determined by the circular velocity at halo virial radius, which is about 2.3 times larger than the one adopted in panel (a). Note that the 
exp(— t v ) model and the full calculation in panel (b) assume different initial Lya line widths, so it is not an apple-to-apple comparison. It 
shows that modifying the initial line width of the exp(— t„) model does not lead to a result mimicking that from the full calculation. 



narrow range, in particular towards higher values of the 
"true" ratio. Even within the narrow range, the ratio 
from the exp(— t v ) model has a large scatter (a factor of 
a few) at a fixed "true" ratio. 

The Lya flux suppression in our exp(— t u ) model ap- 
pears to_j3e_j i nu£h_^ronj^^ work 
(e.g.. lMcQuinn et alJl2007t llliev et al.ll2008| ). Compared 
to our work, previous work assumes a much wider ini- 
tial Lya line width and lower gas temperature (set to be 
10 4 K). The differences in the above two factors largely 
explains the differences in the results. For more detailed 
explanations and discussions, see Appendix [XJ where we 
perform several tests with the exp(— r„) model by vary- 
ing the initial Lya line width and gas temperature. 

Figure [12] compares the correlation between apparent 
to intrinsic luminosity ratio and peak wavelength shift 



from our full radiative transfer model (dotted contours) 
and the above two exp(— r„) models (solid contours). Al- 
though the sign of the correlation is the same for all three 
models, the exp(— r„) model with smaller (larger) intrin- 
sic Lya line width gives a smaller (larger) slope in the 
correlation than the full calculation. 

The comparisons with the exp(— t„) model results in 
this section demonstrate that the exp(— t v ) model can 
provide a qualitative understanding of the results, but 
there is no simple way to modify the results from the 
exp(— Ty) model to match those of the full calculation. 
Full radiative transfer calculation is necessary to obtain 
quantitatively correct results in Lya emission properties 
of LAEs. 

5.4. Lya Luminosity Function of LAEs 



14 



a: O.l 



„ io- 2 



I 10" 3 



10- 





'}"[ 


- log[M/(h- 


■M s )]<10.5 


" log[M/(h- 


■M e )]il0.5 

/ _ 







io- 



10" : 



10" s 0.1 



Fig. 11. — Comparison between the apparcnt-to-intrinsic Ly« lu- 
minosity ratios from our full radiative transfer calculation (full RT) 
and the exp(— t„) model. The thin solid curve shows the median 
ratio from the exp(— t„) model as a function of the ratio from the 
full RT for halos below 1O 1O ' 5 /i _1 M . The two thin dotted curves 
indicate the lower and upper quartiles. The set of thick curves are 
for halos above 10 10 ' 5 /i — Mq. The diagonal dashed line is the line 
of equality. 




log (Lapparent/Ljntrinsic) 

Fig. 12. — Comparison of the correlation between apparent to 
intrinsic luminosity ratio and peak wavelength shift from different 
models. The dotted contours are from our model with full calcu- 
lation of Lya radiative transfer (the same contours as in Fig. 11211 . 
The solid contours on the left are from the simple exp(— t„) model 
with the initial Lyci line width given by halo virial temperature, as 
adopted in our model with full radiative transfer calculation. The 
solid contours on the right are also from the exp(— t„) model, but 
the initial Lya line width is determined by the circular velocity 
at halo virial radius, which is about 2.3 times larger than the one 
from halo virial temperature. For all the three sets of contours, the 
peak wavelength shift is in units of the line width a\ set by halo 
virial temperature. The two contours in each case show the 68% 
and 95% distribution, respectively. 



An important product of surveys of LAEs is the Lya 
LF, one of the most widely studied statistical proper- 
ties of LAEs. The LF can be used to infer the relation 
between LAEs and their host dark matter halos. The 
evolution of LFs around the reioniza tion epoch can probe 
the reionization of the universe fe.g-. IMalhotra fc Rhoadsl 
2004; Ha iman fc Cen1 l2005h As we have shown in S 15.31 
the observed Lya luminosity of LAEs differ from the in- 
trinsic one. Here we study the Lya LF of LAEs from our 
model. 

The star formation prescription adopted in the reion- 
ization simulation leads to a tight correlation between 
star formation rate and halo mass. Therefore the in- 
trinsic Lya luminosity is tied to the halo mass and the 
intrinsic Lya LF largely reflects the halo mass function. 
In Figure 1131 the filled squares show the intrinsic Lya 

LF, $j (^intrinsic 

). The apparent Lya LF, $a(£ a ppare n t), 
is related to the intrinsic one through 



apparent J 



p(L 



apparent | -^intrinsic 



where p(L appa . Tent \L 



(3) 

) is the probability density of 



intrinsic. 

the apparent luminosity at cl given ^intrinsic (Fig. &)■ 
The apparent Lya LF can be directly read off from Fig- 
ure IHJi and is shown as open squares in Figure [T3J In 
terms of luminosity, the apparent Lya LF shifts towards 
the faint end by a factor of 5-20 with respect to the 
intrinsic one and the shift is larger at the faint end. For- 
mally, the intrinsic LF can be inferred from the observed 
one by 
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intrinsic ^apparent) is the probability density of 
the intrinsic luminosity at a given L a pparent (Fig. HJ;) 



where p(L 



However, even if p(Lintrinsic|£ap P arent) is known or as- 
sumed a priori, it is not enough to infer the intrinsic LF 
from the observed one. The reason is that at a given 
intrinsic luminosity, the distribution of apparent lumi- 
nosity can have a tail to low value (Fig. [8j) . So one needs 
to have observations of the apparently faint LAEs to re- 
cover the full information of the intrinsic luminosity dis- 
tribution, or one has to rely on the extrapolation of the 
observed LF to the faint end. 

For comparison, the crosses in Figure [T31 show the Lya 
LF from the exp(— t„) model, which adopts the same 
intrinsic Lya line width as the full radiative transfer cal- 
culation. As already demonstrated in § 15.31 the Lya flux 
is more strongly suppressed in the exp(— t„) model. The 
resultant apparent Lya LF looks like the intrinsic one 
shifting toward the faint end by more than two orders of 
magnitude in luminosity. With the same intrinsic Lya 
LF, the apparent Lya LF from the full radiative transfer 
model is significantly higher than that from the exp(— r v ) 
model, a consequence of the frequency and spatial diffu- 
sion of Lya photons from scatterings. 

6. IMPLICATIONS FOR OBSERVED PROPERTIES OF LAES 

Our model of LAEs predicts the relation between ob- 
served Lya emission and the intrinsic one. With simple 
assumptions about the intrinsic properties of LAEs, we 
are able to predict an array of observational properties of 
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Fig. 13. — Lycv luminosity function of z ~ 5.7 LAEs. The 
intrinsic luminosity function in our model is represented by filled 
squares. The apparent luminosity function from a full radiative 
transfer calculation is plotted as open squares. For comparison, 
the crosses show the apparent luminosity function with a simple 
exp( — Ty ) radiative transfer treatment. Poisson errors are plotted. 
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Fig. 14. — Observed Lyo LF of z ~ 5.7 LAEs. Solid circles 
connected by solid lines are the o bserved LF for z ~ 5.7 LAEs in the 
1 deg 2 SXDS (Ouchi et al. 2008). Open circles connected by dotted 
lines are the observed one from t he 0.2 deg 2 Subaru De ep Field 
(a different sky area from SXDS; IShimagaku et al.1 12006ft . Open 
squares are the apparent LF from our model with a full radiative 
transfer calculation, and diamonds are the same but with the Lya 
luminosity boosted by a factor of 5 (see text for details). 

LAEs. In this section, we compare our model predictions 
to observations for z ~ 5.7 LAEs and attempt to under- 
stand various observed properties of LAEs. We focus 
here on the Lya LF, the UV LF, and the Lya equivalent 
width (EW) distribution of LAEs. 

6.1. Lya Luminosity Function 

In our model, the apparent Lya luminosity corre- 
sponds to the observed one. The open squares in Fig- 
ure [14] (the same as in Fig. [13J) show the predicted 



Lya LF of LAEs from the model. Filled circles are 
the me asurement based o n 401 LAEs from the 1 deg 2 
SXDS (|Ouchi et all 12008ft . while open circles show that 
based on 89 LAEs in a 0.2 deg 2 Subaru Deep Field 
(jShimasaku et all 12006ft that covers a different sky area 
from SXDS. 

Compared with the observed Lya LF, the apparent LF 
from our model appears to be one order of magnitude 
lower in normalization. This seems to put into question 
our model. However, the normalization is not necessarily 
a good indicator of the importance of the difference. A 
more sensible way is to characterize the difference by the 
shift in luminosity scale — with respect to the observed 
Lya LF, the apparent LF from our model is shifted by 
a factor of 3-6 to the left (toward the low luminosity 
end) . If we simply increase the apparent Lya luminosity 
by a factor of 5, the apparent LF (shown as diamonds in 
Fig. [T4|) shows a better match to the observed one. The 
factor of 3-6 discrepancy between the model and the data 
implies that we are missing some physics in our model- 
ing. For comparison, the apparent Lya LF based on the 
simple exp(— r„) model (crosses in Fig. HUl) corresponds 
to about two orders of magnitude shift in luminosity with 
respect to the observed one. Again, the frequency and 
spatial diffusion from the realistic radiative transfer in 
our model enhances the probability of being detected for 
Lya photons. 

Before discussing possible solutions to the factor of 3-6 
shift in luminosity, we justify that changing the normal- 
ization of Lya LF cannot be the solution for the dis- 
crepancy. The normalization can only be changed by 
changing the amplitude of the halo mass function. How- 
ever, we do not expect a large uncertainty in the latter 
with the current constraints on cosmology. The reion- 
ization simulation that our model is based on adopted 
cosmological parameters in accordance with WMAP five- 
year results. To estimate the cosmology-caused change 
in the halo mass function, we make use of the analytic 
formula given bv lSheth fc Tormenl lfl999ft. Not surpris- 
ingly, the uncertainty in as dominates the amplitude 
change in the halo mass function. Within the 1-a un- 
certai nties of the WMAP five-year cosmological param- 
eters (jDunkley et al.ll2009h . the cosmological parameters 
used in the simulation already put the halo mass function 
amplitude in the high end. Different combinations of cos- 
mological parameters within their 1-a uncertainty ranges 
can only boost the amplitude by ~ 30%. Therefore, 
we conclude that cosmology-caused amplitude change in 
halo mass function does not help much in making our 
model Lya LF match the observed one. 

Lya luminosity depends on the amount of massive 
stars that emit ionizing photons. The conversion from 
star formation rate t o intrinsi c Lya luminosity adopted 
in our model assumes ISalpeterl (Il955ft IMF. For a fixed to- 
tal stellar mass with star mass dist ributed in the ra nge of 
0.1-100M Q , simply changing to a lChabrierl ((2003) IMF, 
which has shallower slopes than the Salpeter IMF below 
a characteristic mass of IMq, can increase the ionizing 
photons by a factor of ~1.6. This assumes no evolu- 
tion in t he stellar IMF with cosmic tim e. Recent stud- 
ies (e.g.. iDavel [20081 : Ivan Dokkuml 12008ft show evidence 
of an evolving IMF and observations of z ~7 dropout 
galaxies also implies an IMF changes towards high red- 
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Fig. 15. — Effect of the initial Lya line width. Plotted here is 
the ratio of the apparent (observed) to intrinsic Lya luminosity of 
z ~ 5.7 LAEs as a function of the initial Lya line width Ci n it. The 
bottom axis marks the line width in terms of the width cr v ir set by 
halo virial temperature, while the top axis in units of the circular 
velocity v c at halo virial radius. The solid and dotted curves are 
the median and quartiles of the luminosity ratio distribution. This 
test is done with a small simulation (box size of 25fe — x Mpc on 
a side). The distribution is calculated from all sources with halo 
mass above 5 X 10 9 /i — 1 Mq. See text. 

shift (e.g.. lOuchi et al1 l2009). The form of evolution can 
be thought as an increasing characteristic mass in the 
Chabrier IMF with increasing redshift. By adopti ng ei- 
ther fo rm o f the IMF at high re dshift proposed by iDavel 
(|200l and Ivan DokTulr] ((2008), we find an increase in 
the ionizing photons by a factor of ~5.5 with respect 
to the Salpeter IMF. Therefore, for the instantaneous 
star formation rate, the Lya luminosity could have the 
amount of enhancement we need for the model predic- 
tion to match the observation. Obviously, the situation 
is not as simple as this and there are additional factors, 
such as the dependence on star formation history and 
metallicity. Overall, a possible factor of 3-6 in under- 
estimating the intrinsic Lya luminosities in our present 
model may be physically accounted for and that would 
put our model in accord with observations. 

In addition, a possible increase in the intrinsic Lya 
width may also increase the apparent luminosity. In 
our model, the intrinsic Lya line profile is assumed 
to be Gaussian with width determined by the halo 
virial temperature. The virial temperature is d efined 
as T vir = GM hf im H /(3kR vil ) (|Trac fc Ceri 120071) . re- 
sulting in a line width of (/i/S) 1 / 2 ^ ~ 0.44w c . Here 
fi ~ 0.59 is the mean molecular weight and v c ~ 
155[M/ t /(10 11 /i" 1 Af Q )] 1 / 3 kms _1 the circular velocity at 
the virial radius. Such an assumption on the line width 
is conservative. Many processes can alter the intrinsic 
line profile. Disk rotation can change the line profile and 
broaden it. Galaxy merging and galactic wind from star 
formation can also substantially broaden the line. Many 
authors have adopted line width ~ v c or larger in com- 
puting Lya transmission with the exp(— t v ) model (e.g. , 
Santos! 120041: iDiikstra et~al1l2007t iMcQuinn et~aLll2007t 
Iliev et al. We have tested the effect of the initial 



line width by performing Lya radiative transfer for LAEs 
in a small simulation of box size of 25/i~ 1 Mpc as well as 
for a subset of LAEs in the 100/i _1 Mpc box simulation. 
Figure [15] shows the apparent to intrinsic Lya luminos- 
ity ratio as a function of the initial Lya line width from 
the simulation with the 25/i _1 Mpc box, inferred from all 
sources with halo mass above 5 x lO 9 ft. _1 M0. Changing 
the line width we adopt to v c can lead to a factor of ^5 
increase in the median luminosity ratio. Therefore, the 
effect of initial Lya line width can potentially shift the 
apparent Lya LF by a factor of a few towards the bright 
end. 

Moreover, the star formation prescription adopted in 
the simulation has a number of assumptions and uncer- 
tainties because of the limitation in our understanding of 
the baryon-related processes. A higher SFR than that in 
the simulation would lead to an increase in the intrinsic 
Lya luminosity, and hence an increase in the apparent 
Lya luminosity. 

We see that a combination of changes in stellar IMF, 
intrinsic line width, and SFR can solve the problem of a 
factor of a few mismatch between the apparent Lya LF 
of our model and the observed one. Changing the IMF 
and SFR would also change the reionization process and 
the neutral fraction of IGM, although they are degener- 
ate with uncertain escape fraction of ionizing photons. 
The prediction of Lya luminosity is coupled with the 
evolution of the gas ionization state. A self-consistent 
calculation is possible, but it is out of the scope of this 
paper. An accurate estimation of the initial Lya profile 
and width would require detailed calculation of the Lya 
transfer through realistic ISM of high redshift galaxies, 
which is little constrained and difficult to compute from 
first principles presently. We conclude that, owing to 
model uncertainties, the discrepancy between model LF 
and observed LF is not as serious a problem as it ap- 
pears to be and may be indications of some interesting 
physics that is not considered in our current calculation 
or the need to have more accurate prescriptions of some 
processes. 

6.2. UV Luminosity Function 

As mentioned in £|5.3[ at a fixed intrinsic Lya lumi- 
nosity of LAEs, the apparent (observed) Lya luminosity 
has a broad distribution, and vice versa. Since the in- 
trinsic Lya luminosity is directly correlated with the UV 
luminosity, our results mean that the UV LF of the ob- 
served LAEs (with apparent Lya luminosity above cer- 
tain threshold) must differ from the intrinsic one. In 
what follows, we show the differences caused by the Lya 
selection and compare the model UV LF of LAEs with 
those from observations. 

We convert the SFR in halos to the UV luminosity Ltjv 
(at 1500 A) through 

Luv - 8 x 10 27 [SFR/(Af Q yr- 1 )] ergs" 1 Hz" 1 , (5) 

which assumes Sal peter IMF and solar metallicity 
(jMadau et al.l Il998| ). Following observers, we ex- 
press the UV luminosity in AB magnitude, Mtjv = 
-2.51og[L uv /(47rdg)] - 48.60 with d = lOpc. 

The thick solid curve in the top- left panel of Figure [T5] 
is the UV LF for all sources (galaxies) in our model. 
Because of the tight correlation between SFR and halo 
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Fig. 16. — UV LFs of LAEs. Left panels are for z ~ 5.7 LAEs. The thick solid curve is the UV LF for all galaxies in our model. The 
filled squares connected by the thin solid curve are the UV LF of model LAEs with the apparent Lyo luminosity above a threshold. The 
threshold is chosen so that the number den s ity m atches tha t of SXDS z = 5.7 LAEs (see text) . The observed UV LF s of LAEs are plotted 
as open symbols, obtained by[Ouchi ct al. (2008) (circles), Shimasaku ct al. (2006) (triangles') . andlHu et aLl ||2006 | ) (stars), respectively. 
The LF of i-dropout galaxies at z ~ 6 is shown as crosses (Bouwcns et al. 2006) and asterisk (Shimasaku et al. 2006]). In the bottom left 
panel, an effective UV extinction of Ajjv e ff = 0-3 is applied to the two curves from the model (see text). The two right panels are the 
observed UV LF for LAEs at 2 ~ 3 and z'~ 4 in SXDS, taken from lOuchi et alj (I2008T) . 



mass, the curve is basically a transformation of the halo 
mass function with a constant mass-to-light ratio. This 
UV LF is from all the galaxies in our model, regardless of 
the apparent Lya luminosity. To be detected as LAEs, 
the apparent Lya luminosity should be high, which im- 
poses a selection function onto the full UV LF. The UV 
LF of LAEs is from sources with apparent Lya luminos- 
ity above a threshold, which is set in our model such that 
the number density of the selected LAEs (that can be 
observed) matches that of the z = 5.7 LAEs in SXDS 
(about 4.0 x 10~ 4 Mpc~ 3 for our adopted cosmology). 
The filled squares connected by thin solid curve shows 
the UV LF for the LAEs selected this way. For a cut in 
apparent Lya luminosity, sources with a higher intrinsic 
Lya luminosity (hence a higher UV luminosity) have a 



higher possibility to be selected (see Fig|5]). As a result, 
the UV LF of LAEs is close to the full UV LF at the high 
luminosity end. However, in lower mass halos (sources 
with lower UV luminosity), fewer sources can have ap- 
parent Lya luminosity high enough to be detected. The 
UV LF of LAEs becomes lower than the full UV LF. 
The ratio of the UV LFs of LAEs and all galaxies as a 
function of UV luminosity is nothing more than a reflec- 
tion of the distribution of apparent Lya luminosity as a 
function of halo mass (FigJSJ) . As a consequence of this 
simple effect, the UV LF of LAEs becomes flattened as 
the UV luminosity decreases and drops rapidly towards 
the low luminosity end (Fig fTTJ)) . 

The predicted features in the UV LF of LAEs are 
indeed seen in observations. T he cyan points are ob- 
served UV LFs of z = 5.7 LAEs (jShimasaku et alj|2005 
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IHu et al.|[200llOuchi et al.ll2008l ). The LF becomes flat 
for Muv > —20.5. The z — 5.7 observations are not 
deep enough to show the predicted drop of the UV LF 
at the faint end, but the drop can be clearly seen in the 
z ~ 3 and z ~ 4 data, as shown in the righ t panels of 
FigHS (also see Fig.22 of lOuchi et all 120081 ) ■ Without 
any adjustment, the UV LF of LAEs from our model is 
in quite a reasonable agreement with observations. The 
observed flattening of the LF toward lower luminosity is 
well explained by our model. The agreement improves 
by adding an effective UV extinction of ^4uv,cff = 0.3 
to our model curve (lower left panel of FigJTBJ. As im- 
plied in § 16.11 the assumed IMF and the adopted SFR 
in the simulation may not be accurate, and we may need 
a higher UV luminosity. Therefore, the effective UV ex- 
tinction here should be understood as a combination of 
t he model uncertainty an d dust extinction. 

iKobavashi et alJ (|2010( ) present an LAE model with 
Lya escape fraction and UV extinction being functions 
of metal column density and starforming and outflow 
phases of galaxies, which roughly repr oduce the observed 
UV LF. The semi-analytic model of (jSamui et al.ll2009f ) 
with constant Lya escape fraction and UV extinction 
fails to reproduce the turnover of the UV LF towards 
low luminosity end. By contrast, there is no Lya escape 
fraction parameter and mass dependent UV extinction in 
our model, and the radiative transfer is the single factor 
responsible to convert the intrinsic Lya emission to the 
observed one. In other words, the Lya escape fraction, 
defined as the ratio of apparent to intrinsic Lya lumi- 
nosity is an output of the model. It is encouraging that 
our model, by accounting for simple physics, is able to 
reproduce the features in the observed UV LF. This is 
an independent output of our model and clearly lends 
credence to our model. 

The sources that are not detectable as LAEs because of 
a low apparent Lya luminosity can be de tected as galax- 
ies th rough the dropout technique (e.g., iBouwens et all 
I2003D . In Figure [TH the observed UV LF of the z ~ 6 i- 



dropo ut galaxies is shown as red crosses (|Bouwens et alJ 
120061 ) and asterisk (jShimasaku et al.ll2006| ). At the faint 
end, it agrees with the full UV LF from our model. It falls 
slightly steeper than our model curve at the bright end. 
Note that the i-drop out technique can miss LAEs with 
strong Lya emission (jOuchi et al.ll2008t ). The face values 
of UV LFs from our model and the data suggest that the 
sum of UV LFs of LAEs and dropout galaxies makes the 
full UV LF. Obviously, the sum should not double-count 
those LAEs that are also detected as dropout galaxies. 
Observations show that abou t 30-50% of dropout galax- 
ies are detected as LAEs Ce.g..lDow-Hygelund et al.ll2007l : 
IStanwav et al.l 120071 lOuchi et alJ 120081 ) (our model im- 
plies that this fraction can depend on UV luminosity). 
Roughly accounting for this fraction, the nonduplicated 
sum of the observed UV LFs of LAEs and i-dropout 
galaxies seems to be in a reasonable agreement with the 
full UV LF from the model. This is yet another indepen- 
dent output of our model that agrees with observations. 

We see that, as a consequence of Lya radiative trans- 
fer, LAEs are sources that have a strong Lya selection 
imposed. This selection effect nicely explains the shape 
of the observed UV LF of LAEs and that of the dropout 
galaxies. 
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Fig. 17. — Distribution of rest-frame Lya equivalent width 
(EWo) from the model. Different line types correspond to different 
UV magnitude. 

6.3. Distribution of Lya Equivalent Width 

As shown in the above two subsections, with simple 
assumptions to account for the model uncertainty, the 
model is able to reproduce the Lya LF and UV LF of 
LAEs. We can go beyond the two LFs to study the joint 
distribution of Lya and UV luminosities, which can be 
casted as the distribution of rest-frame Lya EW as a 
function of UV luminosity. This distribution is not lim- 
ited to LAEs and it can include that from the dropout 
galaxies. 

Observationally, it is found that z ~ 6 galaxies seem 
to show a d eficit of large EW values for UV lumi- 
nous objects flAndo et al.ll20"ol iShimasaku et al.l [20061: 



IStanwav et all 120071: lOuchi et al.ll2008[ ) (see the data 
points in Fig. [T8|) . The threshold UV luminosity for th e 
deficiency is M uv = -21.5 to -21.0 (jAndo et alJl2006f ). 
Similar trend is seen for z ~3-5 L AEs as well (e.g., 
Ouchi et all 120081: iShiova et al J 12009). The trend is also 
reported for dropout galaxies at different redshifts (e.g 



Shaplev et alj|2003t lAndo et al.ll2007t|Kaiino et alj|2009t: 
Pentericci et al.ll2009HVanzella et al.ll2009D . lAndo et alJ 
(2006) and other authors invoked the differences in dust 
extinction, amount of internal and surrounding neutral 
hydrogen gas, age of stellar population, and/or gas kine- 
matics between UV faint and luminous galaxies as pos- 
sible c auses of the trend seen in Lya EW and UV lumi- 
nosity. iMao et al. (2007) present a model of high redshift 
gala xies including chemical evolution and dust attenua- 
tion. Kobayashi et al. (2010) also present a semi-analytic 
model of LAEs, in which Lya and UV are attenuated dif- 
ferently by clumpy dust distribution. Both models seem 
to explain the deficiency of high Lya EW in luminous 
galaxies essentially by means of halo mass dependent 
dust content. On the other hand, based upon a pure 
statist ical analysis with z ~ 3 galaxies, iNilsson et alJ 
conclude that there is no dependence of Lya EW 



on UV luminosity for LAEs and Lyman break galaxies 
(LBGs). They interpret the lack of large Lya EW, UV 
bright galaxies as an observational effect of small survey 
volumes. 
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Fig. 18. — Rest-frame Lye? equivalent width (EWo) as a function of UV lumino sity. In both panel s, the thick solid curve corresponds 
to the threshold of the observed (apparent) Ly q luminosity f or the z ~ 5.7 LAEs in lOuchi et al.l 1120081 ), and the thick dashed vertical line 
shows the 3— a limit for the UV photometry in Ouchi ct al. (2008). Dotted curves in both panels show the distribution of objects in the 
plane of EWo and UV luminosity from our model (see text): in the left panel, the contours denote the probability density of objects in the 
plane (thicker contours for higher den sities); in the right panel, we plot the median (thin solid curve) and quartiles (dotted curves) of EWo 
for the LA Es that can be dete cted in Ouchi et al. (2008), i.e., for sources above the thick solid curve. The data points in both pane ls are 
taken from Ouchi ct al. (2008): the filled circles and squares are for z ~ 5.7 LAEs in Ouchi et al. (2008) and Shimasaku ct al. (200( f) from 
the Su baru fields; the stars and diamonds represent z ~ 6 dropout g alaxies in Stanway ct al. (2007) and those compiled bv lAndo et al.l 
(2006); the crosses denote z ~ 5 dropout galaxies in Ando ct al. (2006). 



For the model presented in this paper, the UV lumi- 
nosity is directly related to halo mass and the appar- 
ent Lya luminosity is determined by radiative transfer, 
which depends on the environment of galaxies (§EJ). It 
is interesting to see to what extent the observed relation 
between Lya EW and UV luminosity can be explained 
by our model. Following § 16.11 and § 16.21 we scale the 
Lya luminosity by a factor of 5 and apply an effective 
UV extinction of 0.3 for the UV luminosity for each LAE 
source. 

In our model, the intrinsic Lya EW is a constant for 
all sources, since both the intrinsic Lya luminosity and 
the UV luminosity are proportional to the SFR. Radia- 
tive transfer, however, gives rise to a broad distribution 
of the apparent (observed) Lya EW. The distribution of 
the apparent EW is similar to that in Figure |5f>, if the 
horizontal axis is relabeled. Figure [8f> shows the distri- 
bution in logarithmic space, while EW distribution in 
linear space matches more closely with what can be in- 
ferred from observations. In Figure |T7l we show the dis- 
tribution of rest-frame Lya EW in linear space from our 
model. At a given UV luminosity, the distribution func- 
tion of apparent Lya EW is a decreasing function of EW. 
In the UV luminosity range considered here, the distri- 
bution function drops faster for sources with lower UV 
luminosity. 

In the left panel of Figure [TSJ the cyan contours show 
the probability density distribution of objects in the 
plane of Lya EW and UV luminosity from the model 
(thicker contours for higher densities) . The apparent Lya 
EW distribution at a fixed UV luminosity roughly follows 
an exponential distribution from our model, with fewer 
sources having larger EW. If the EW distribution did not 
vary with UV luminosity, the contour of equal probabil- 
ity density in the left panel of Figure [T8l would appear to 
tilt along the direction of low EW and UV bright to high 



EW and UV faint because the number density of objects 
drops fast with UV luminosity. In our model, the ap- 
parent EW distribution has a weak dependence on UV 
luminosity, but it only leads to a slight change in tilt di- 
rection of the contour and cannot mask the effect caused 
by the decreasing number density toward high UV lumi- 
nosity. The cyan contours in the left panel of Figure [15] 
clearly show that sources at the corner of large appar- 
ent EW and high UV luminosity have a low probability 
density. The low probability is a consequence of the com- 
bination of two facts: that UV LF drops steeply towards 
high luminosity and that the distribution of the apparent 
Lya EW at fixed UV luminosity is a decreasing function 
of EW. The result suggests that a large survey volume is 
needed to discover large EW, UV bright sources , whic h 
is consistent with the conclusion in lNilsson et al. (2009). 
We note that we do not assume any particular form of the 
EW distribution at a given UV luminosity. In our model, 
we have a single value of intrinsic EW for all sources, 
since both intrinsic Lya luminosity and UV luminosity 
are proportional to SFR and we apply the same scaling in 
either luminosity for all sources. The distribution of the 
apparent EW simply results from the radiative transfer 
effect. 

Because of the flux limit, luminosities of observed 
LAEs are above a threshold. The red solid curve in 
either panel of Figure [TH is the corresponding thresh- 
old for EW as a function of UV luminosity for z = 5.7 
LAEs in SXDS. To make a further comparison between 
the model and the observed LAEs, we show in the right 
panel of Figure [TS] the median and quartiles (solid and 
dotted cyan curves) of the EW distribution as a function 
of UV luminosity from the model, for LAEs correspond- 
ing to those observed in the Subaru fields (filled circles 
and squares). The model prediction largely follows the 
observational trend. Again the lack of large EW, UV 
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Fig. 19. — Distribution of Lya rest-frame equivalent width 
(EWo) for z ~ 5.7 LAEs. The dashed curve and the dotted his- 
togram are two estimates of the distr ibution for all the photomet- 
rically selected z ~ 5.7 LAEs in IQuchi etaLl [120081) . The dashed 
curve is estimated from a maximum likelihood method by fully ac- 
counting for the probability distribution of the measured EWo for 
each LAE and the dotted histogram is o btained by si mply adopt- 
ing the best measurement of EWo ( scc Ouchi ct al. 2008). The 
solid histogram is from our model for the corresponding popula- 
ti on of LAEs, i . e., for LAEs that have the same number density as 
in lOuchi etU] 1I2OO8I ). 

bright LAEs is evident in the model. Although the small 
number of observed LAEs prevents a comparison of the 
model and observed EW distribution as a function of UV 
luminosity, we can make a comparison for the overall EW 
distribution of the observed LAEs. The dashed curve 
and dotted histogram in Figure [19] are two estimates of 
the EW di stribution for all th e LAEs detected in SXDS 
inferred bv lOuchi et all pOOl . The dashed curve is ob- 
tained with a maximum likelihood method by accounting 
for the full probability distribution of the measured EW 
for each LAE. The dotted histogram is obtained by sim- 
ply counting the number of LAEs in each EW bin based 
on the measured values of EW, i.e., no u ncertainty in 
the m easured EW is assumed. According to lOuchi et al.l 
(2008), the two estimates likely bracket the true distri- 
bution. The corresponding EW distribution from our 
model is shown as the solid histogram, which appears 
to be in good agreement with the observation estimates. 
Interestingly, it is more closely resemble the one from the 
maximum likelihood method. 

Without appealing to differences in UV faint and UV 
bright sources, such as the amount of dust and age of stel- 
lar population, the observed deficit of bright UV galaxies 
with large Lya EW is reproduced by the model. It is a 
natural consequence of the fact that UV LF drops toward 
high luminosity and the distribution of apparent EW at 
fixed UV luminosity is largely a decreasing function of 
EW. The model also reproduces the observed EW distri- 
bution. We emphasize that at fixed UV luminosity, the 
distribution of apparent EW in the model is completely 
determined by Lya radiative transfer effect, not by any 
other mechanisms (e.g., different UV extinction or stel- 
lar age). Our model suggests that dependence of dust or 
stellar age on halo mass (or UV luminosity) , if there is 



any, does not play a dominant role in the observed Lya 
EW distribution of high-redshift galaxies and that Lya 
radiative transfer is the main mechanism in determining 
the observational properties. 

7. IMPORTANT PHYSICAL FACTORS IN DETERMINING 
THE OBSERVABILITY OF LAES 

In previous sections, we study statistical properties of 
Lya emission of LAEs in our model. At a fixed intrinsic 
Lya luminosity, the apparent luminosity and peak wave- 
length shift have broad distributions. The cause of the 
distribution must be related to the underlying distribu- 
tion and kinematics of gas, hence the matter distribution, 
around LAE sources, since Lya radiative transfer is sensi- 
tive to the density and velocity fields. In this section, we 
study the correlations between the apparent Lya prop- 
erties and the environment of matter around sources. By 
revealing the physical causes, such correlations would aid 
our understanding of the observability and the statistical 
properties of LAEs. 

To identify the key factors in shaping the observability 
of LAEs, we first choose LAE host halos in a narrow mass 
bin (~ 8 x 10 10 /i _1 M Q ) and stack the neutral gas density, 
temperature, and peculiar velocity profiles along the line 
of sight centered on these halos. As shown by the dotted 
curve in the middle-left panel of Figure [201 the stacked 
density profile around a halo appears to peak at the halo 
center and to be symmetric around the center. The cen- 
tral density is about a factor of 2000 higher than the 
mean cosmic density. We then stack only sources that are 
strongly and weakly suppressed in Lya luminosity (e.g., 
the upper and lower quartiles of Lapparcnt/^intrinsic), re- 
spectively, into two subsets. The stacked density profiles 
of strongly and weakly suppressed sources (thin and thick 
solid curves in the middle-left panel of Fig. [20]) appear 
to have a relative offset in amplitude and are asymmet- 
ric in opposite directions. The trend becomes clear with 
the total mean profile subtracted (middle-right panel of 
Fig. l20[) . The temperature profiles of strongly and weakly 
suppressed Lya sources also show systematic but small 
differences (top panels of Fig. [20)) . As temperature is 
largely linked to density inside halos, we tend to incor- 
porate its effect into an overall density effect. The differ- 
ence in amplitude and asymmetry of density profiles be- 
tween strongly and weakly suppressed sources indicates 
that density and density gradient along the line of sight 
contribute to the observability of LAEs. 

The overall stacked peculiar velocity profile around a 
halo (bottom- left panel of Fig. |2"U)) . with the velocity cen- 
tered on the halo velocity, shows clear signatures of infall 
region: an inner Hubble-like contraction region near the 
center and an outer region with infall velocity decreasing 
outward. Similar to the density profile, we find differ- 
ence in the slope of the velocity profile in the outer in- 
fall region between strongly and weakly suppressed Lya 
sources (bottom- right panel of Fig. |2"0|) . Therefore, the 
peculiar velocity gradient along the line of sight is a fac- 
tor in determining the observability of LAEs. We check 
the dependence on the line-of-sight velocity of halos and 
find that halo velocity also contributes to the LAE ob- 
servability. 

The above exercises provide us with the initial evi- 
dence on what physical variables affect the observability 
of LAEs. Since on scales larger than halo size, gas den- 
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Fig. 20. — Mean gas distribution profiles along the line of sight around halos of mass 8 X W 10 h~ 1 Mq. In the left panels, from top 
to bottom, are temperature, neutral hydrogen density, and peculiar velocity profiles. Note that the halo velocity is subtracted from the 
peculiar velocity. The dotted curve in each panel is the mean profile for all the halos in the narrow mass bin around 8 X 10 10 /i _1 Mq. The 
thin and thick solid curves are those in the lower and upper quartiles of the apparent to intrinsic Ly« luminosity ratio. In the corresponding 
panel on the right, we show the profile with respect to the mean profile, i.e., the mean profile for all halos in the mass bin is subtracted. 



sity and velocity largely follow those of the underlying 
dark matter, we identify the matter density and peculiar 
velocity and their gradients along the line of sight as the 
major factors in shaping the LAE observability. We pro- 
ceed to study the correlation between these quantities 
and the suppression in Lya luminosity for all the sources 
in our model. We arc interested in the environment den- 
sity and velocity around sources, not those inside halos, 
but they cannot be obtained directly at the positions 
of halos from the outputs of the reionization simulation. 
To eliminate the influence of density and velocity profiles 
inside halos, we smooth the density field with a 3D top- 
hat filter of radius of 2/i _1 Mpc (comoving). Conclusions 
reached below are largely immune to possible uncertain- 
ties related to limited resolution of the hydro simulations 
we use (^65kpc comoving), which nonetheless is much 
smaller than 2/i _1 Mpc comoving. 



With the smoothed overdensity field 5, we solve the 
linear peculiar velocity from the continuity equation 



(5+iv-v 
a 



0, 



(6) 



where 5 can be written as fH(a)S with / = din 5 /din a 
the derivative of the growth factor and H(a) the Hubble 
parameter at the time when the scale factor is a. The 
peculiar velocity field is obtained from the Fourier trans- 
form 



(7) 



where 5k is the Fourier transform of the smoothed over- 
density field. The velocity gradient and density gradient 
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Fig. 21. — Dependence of Lya flux suppression of LAEs on density and peculiar velocity. The suppression is characterized by the ratio 
of the apparent (observed) and intrinsic Lya luminosity Lapparent/£mtrmsic- Panel (a): dependence on the smoothed overdensity field 
at the source position. The overdensity field is smoothed with a 3D top-hat filter of radius 2h~ x Mpc (comoving), which is chosen to be 
larger than the size of the infall region around halos. Panel (b): dependence on the density gradient along the Z direction. The derivative 
is with respect to comoving coordinate. Panel (c): dependence on the host halo velocity. Panel (d): dependence on the linear peculiar 
velocity gradient along the Z-direction. The linear peculiar velocity is obtained from the smoothed overdensity field based on the continuity 
equation (see text for detail). The velocity gradient is put in units of the Hubble parameter. In each panel, the solid curve is the median 
ratio and the two dotted curves delineate the upper and the lower quartiles. Note that the observation is along the — Z direction, which 
matters for interpreting the results in of panels (b) and (c). See the text. 



along the Z direction are 

k 

and 

k 

The spatial derivatives in all the above equations are with 
respect to comoving coordinates. 

Figure [2U shows the correlations of the apparent to 
intrinsic luminosity ratio with the environment overden- 
sity, velocity and their gradient along the line of sight for 



all the LAEs in our model. The trends of correlations 
at a fixed halo mass look similar, with amplitudes and 
slopes of curves slowly evolving with mass. These corre- 
lations reflect different aspects of Lya radiative transfer. 
In what follows, we interpret them in turn. 

Figure [2lb , shows that Lya emission from LAEs is more 
strongly suppressed in denser region. This seems easy to 
understand — higher density means high optical depth 
for Lya emission. At the high end of the density, high 
optical depth also leads to large Lya frequency diffusion 
(and more spatial diffusion), which appears to be able 
to compensate for otherwise reduced Lya transmission 
due to high optical depth. As a result, the curve of the 
overall luminosity suppression factor flattens at very high 
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density. In details, the above understanding is not com- 
plete. Through the continuity equation (eq.[6]), density 
is anti-correlated with velocity gradient. As shown be- 
low (Fig. l2"TTi). velocity gradient plays a dominant role 
in determining the suppression factor. As will be elab- 
orated in Paper II, the density dependence seen here is 
largely driven by its anti-correlation with the velocity 
gradients, while velocity gradients along different direc- 
tions contributes a lot to the anisotropic Lya emission 
distribution. 

Figure l2~Tb shows that Lya sources are easier to be 
transmitted if the local density has a negative gradient 
along the line of sight, i.e., density decreases towards the 
observer. Note that the observation direction (line of 
sight) for Lya emission in our radiative transfer model- 
ing is set to be along the — Z direction, therefore neg- 
ative gradient means that sources are located on the 
near side of overdense regions or on the far side of un- 
derdense regions. With respect to sources located on 
the far side of overdense regions or on the near side of 
underdense regions, Lya photons encounter a decreas- 
ing density profile, thus low optical depth, along their 
way to the observer. Although these photons also travel 
across other overdense and underdense regions, the red- 
shifting of these photons caused by Hubble expansion 
would make these intervening regions largely transpar- 
ent to them. That is, only the immediate environments 
of the sources affect their observability. 

Figure [2Tb show that, with respect to the Hubble flow, 
LAEs moving away from the observer have a lower sup- 
pression in Lya luminosity than those moving towards 
the observer. Note that sources with negative velocities 
are the ones moving away from the observer, since we ob- 
serve along the — Z direction. The IGM on scales larger 
than LAE host halos can detach from the motion of halos 
to some degree, and experiencing Hubble expansion. The 
motion of the halo with respect to the surrounding IGM 
thus introduces a dipole in the Lya optical depth around 
the source. For sources moving away from the observer, 
Lya photons acquire an additional redshift from the halo 
motion, and hence the direction towards the observer is 
the one that has the lower optical depth and photons 
preferentially "leak" towards that direction. 

Figure |2"TH shows that the line-of-sight gradient of the 
line-of-sight velocity has the largest effect in shaping the 
LAE observability. We see that sources located at places 
that have larger line-of-sight gradient in peculiar veloc- 
ity are easier to be observed. A local velocity gradient 
effectively changes the local Hubble expansion rate. A 
positive gradient increases the local expansion rate. A 
faster expansion makes Lya photons on the red side of 
line center much easier to escape. It also makes Lya 
photons on the blue side of line center travel a shorter 
distance to redshift to the line center and to be scattered 
in the IGM, and therefore the escaping photons after fre- 
quency diffusion are more centrally distributed, leading 
to higher surface brightness. Both of these effects cause 
the transmission of Lya photons that is enhanced for 
sources with positive local velocity gradient. 

We see that the major factors in determining LAE ob- 
servability all have clear physical origins. Quantities like 
the density gradient, velocity, and velocity gradient are 
statistically interrelated. For example, in a statistical 
sense, sources on the near side of an overdense region 



(negative density gradient) usually moving away from us 
(negative velocity). However, on a source by source ba- 
sis, because of the randomness of the density and velocity 
field, this is not always true. In other words, there are 
large scatters among the correlations of these quantities. 
Since these different quantities have different physical ef- 
fects on Lya transmission, the overall Lya transmission 
or luminosity suppression effect should be a supposition 
of all of them. 

The dependence of Lya radiative transfer on large scale 
density and peculiar velocity fields imposes a strong se- 
lection effect on observations of LAEs. The selection 
leads to new features in the clustering of LAEs, which 
are investigated in Paper II. 

Some of the environment factors identified here are ex- 
pected to be foun d in the exp(— t„) model. For example, 
llliev et ahl (|2008[) compare the results from the exp(— r„) 
model with and without the peculiar velocity field turn- 
ing off. They find that both the gas infall around halos 
and source peculiar velocity are important in determining 
the suppression factors of Lya luminosity and in shaping 
the observed Lya luminosity, especially at the luminous 
end (see their Figures 21, 23 and 24). This is qualita- 
tively in parallel with our finding that source peculiar ve- 
locity and peculiar velocity gradi ent of surrounding mat- 
ter are important factors. While iMcQuinn et al.l (|2007[ ) 
consider the cases of varying the global neutral fraction, 
which shows the effect of density on the suppression of 
Lya luminosity, it is clearly different from what we in- 
tend to do here. Nevertheless, we expect that the density 
and density gradient effects also show up in the exp(— r„) 
model, given that they contribute to the line-of-sight op- 
tical depth. Although in most cases the simple exp(— r„) 
model can provide qualitative understanding of Lya ra- 
diative transfer results, we do not expect it to capture all 
the physics. In the exp(— t v ) model, the radiative trans- 
fer is completely determined by the line-of-sight optical 
depth. In the full calculation, this is not the case - the 
scatterings of Lya photons enable them to probe the op- 
tical depth in all directions and the line-of-sight outcome 
depends on not only the line-of-sight optical depth but 
also those in other directions (see more details in Paper 
II). 

8. SUMMARY AND DISCUSSION 
8.1. Summary of Main Results 

We perform a full Lya r adiative transfer calcula- 
tion w ith a Monte Carlo code (Zheng_& Miralda- Escudel 
2002) to study LAEs in a cosmological volume. The LAE 
sources and the physical properties of neutral hydrogen 
gas are taken from the z ~ 5.7 outputs of a cosmological 
reionization simulation (jTrac et al.l [2008) . which solves 
the coupled evolution of the dark matter, baryons, and 
ionizing radiation in a box of 100/i _1 Mpc (comoving) 
on a side. The large volume of the simulation allows a 
statistical study of z ~ 5.7 LAEs. Radiative transfer 
of Lya photons in the IGM environment around LAEs, 
which leads to both frequency and spatial diffusion of 
Lya photons, turns out to play a crucial role in deter- 
mining the observability of LAEs and in understanding 
the observed properties of LAEs. 

Although the radiative transfer calculation is compu- 
tationally costly, the LAE model we present in this paper 
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is rather simple. The UV or intrinsic Lya luminosity is 
assumed to be proportional to the SFR, which is tightly 
coupled to halo mass in the reionization simulation we 
use. That is, we essentially adopt a constant mass to 
light ratio, where mass is halo mass and light is cither 
UV or Lya . All we do is to add the physics of Lya 
radiative transfer into the model to obtain the observed 
properties of LAEs. That is, we introduce the radiative 
transfer of Lya photons in the IGM as the single factor 
responsible for transforming the intrinsic Lya emission 
properties to the observed ones. Our model produces 
IFU-like data cube that covers the extent of the simu- 
lation box, which allows mock observations to be made. 
With the Lya image contracted from this d ata cube, we 
follow typical observational procedures (e.g.. lOuchi et al.1 
2008) to identify LAEs and then extract their Lya spec- 
tra. 

Initially Lya photons are produced inside the star for- 
mation region. Therefore the intrinsic Lya sources are 
expected to be similar in size as the UV sour ces, which 
are compact (< 1 kpc; iTaniguchi et al.l [20091. We find 
that an intrinsically point-like Lya source becomes ex- 
tended as a consequence of resonant scatterings of Lya 
photons (spatial diffusion) . The scatterings of Lya pho- 
tons do not destroy them and all Lya photons escape in 
the end. However, observationally, only the central part 
of the extended source can be detected as a consequence 
of the limit set by the surface brightness threshold. The 
scatterings of Lya photons also cause the frequency of 
Lya photons to change (frequency diffusion). The resul- 
tant Lya spectra from the central aperture do not have a 
simple relation to the initial profile, which is assumed to 
be Gaussian in our model. Our results from full radiative 
transfer calculations show a clear difference from a simple 
treatment of Lya radiative transfer, namely the exp(— t„) 
model, widely adopted in previous work, which modifies 
the intrinsic Lya spectrum by multiplying the line-of- 
sight transmission determined by the optical depth at 
each frequency. 

The observed Lya spectrum of an LAE in our model 
shows a clear asymmetry, skewed towards red. Although 
the exp(— t„) model produces the same qualitative fea- 
ture, the predicted line profile, the frequency shift, and 
the total flux are all significantly different from our re- 
sults. While the spectrum of the exp(— r„) model we 
present is essentially the intrinsic one truncated below 
a certain wavel ength (but see Figures 14 and 15 in 
llliev et all 120081 for more complex line shapes, probably 
caused by different assumptions in the exp(— t„) model), 
the observed Lya spectrum in our model can have contri- 
butions from photons with frequency much redder than 
initial photons, a result of the scattering-caused fre- 
quency diffusion. We find that the redward shift of the 
Lya line induced by radiative transfer is usually a few 
times the intrinsic line width, with a distribution that 
peaks at about three times. The asymmetry and shift of 
the Lya line do not indicate the presence of any winds, 
but they arise from the structure of the halo infall and 
Hubble expansion around the sources. If one were to in- 
fer the wind velocity, if there is any, from comparing the 
relative shift in the Lya line and an optically-thin line, 
one has to keep in mind the Lya radiative transfer effect. 
For example, the observationally inferred velocity of the 
receding winds would be overestimated by ~ 100km s _1 



or more if the effect is not taken into account. 

As a consequence of the frequency diffusion and spa- 
tial diffusion, our model predicts a much higher observed 
Lya flux than the exp(— r„) model. At a fixed intrinsic 
Lya luminosity (i.e., fixed host halo mass), the observed 
(apparent) luminosity is broadly distributed. The shift 
in the Lya line peak and the ratio of the apparent to in- 
trinsic Lya luminosity appear to be anti-correlated. The 
distributions of the line peak shift and the ratio of appar- 
ent to intrinsic Lya luminosity, and their correlation, all 
result from the dependence of the Lya radiative transfer 
on the IGM environments around sources. It would be 
interesting and extremely useful if we could make use of 
the full information in the observed Lya properties to 
infer the intrinsic ones, and we reserve such an investi- 
gation for future work. Although our model predicts a 
much higher observed Lya flux than the exp(— t v ) model, 
it still leads to a highly suppressed Lya flux, compared 
with the intrinsic one. The suppression factor depends 
on the assumed line width of the intrinsic Lya spectra. 
For the line width assumed in our model (given by halo 
virial temperature), we find that, with respect to the in- 
trinsic Lya LF of LAEs, the observed (apparent) Lya 
LF shift towards the low luminosity end by roughly one 
order of magnitude in luminosity. For comparison, the 
exp( — t u ) model would shift by two orders of magnitude 
in luminosity. 

We make comparisons between the z ~ 5.7 LAEs in ou r 
model and those observed in SXDS (jOuchi et alj|2008j) . 
The sizes, morphologies, Lya line profiles of the model 
LAEs are remarkably similar to the observed ones. For 
the Lya LF, UV LF, and Lya EW distribution, our 
model can successfully reproduce the observations and 
provide physical explanations for various observed fea- 
tures. 

After an overall adjustment of a factor of ~ 5 in lumi- 
nosity, the Lya LF of model LAEs matches well with ob- 
servation. The adjustment reflects our incomplete knowl- 
edge in the stellar IMF at high redshift, the uncertainty 
in the model SFR, and the lack of information on the 
intrinsic Lya line profile. According to our model, there 
is no one-to-one map between the intrinsic and the ob- 
served Lya luminosity. In other words, there is a large 
scatter in the relation between the apparent and intrin- 
sic luminosities. At a fixed observed luminosity, LAEs 
can differ by one order of magnitude in the intrinsic lu- 
minosity (Fig. |SJ;). This large scatter has to be taken 
into account when interpreting the observed Lya LF and 
linking the observed LAEs to their host halos. 

For the UV LF of observed LAEs, our model prediction 
shows a good agreement with observation. In particular, 
the turnover of the UV LF towards the low luminosity 
end seen in high-z (z ~ 3-6) LAEs are well reproduced. 
The key to interpret the shape of the UV LF is that ob- 
served LAEs are sources with observed (apparent) Lya 
luminosity above certain threshold. The turnover reflects 
that for LAEs with low UV luminosity (or low intrinsic 
Lya luminosity, or low halo mass), the probability for 
the observed Lya luminosity to exceed the observation 
threshold is low, a consequence of the broad distribution 
of apparent Lya luminosity at a given intrinsic Lya lu- 
minosity. The full UV LF for sources in our model (i.e., 
without imposing the observation Lya luminosity thresh- 
old) agrees well with the nonduplicated sum of the ob- 
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served UV LFs of LAEs and i-dropout galaxies at z ~ 6. 

The observed distribution of Lya EW as a function 
of UV luminosity is also reproduced in the model. We 
note that in our model all the sources have the same 
intrinsic EW and the distribution of the observed val- 
ues of EW at fixed UV luminosity is purely caused by 
the environment-dependent radiative transfer effect. At 
a fixed UV luminosity, the distribution of observed (ap- 
parent) EW is a decreasing function toward high values. 
The observat ional trend of lac king UV bright, high EW 
sources (e.g. JAndo et al.ll2006l ) is naturally explained by 
our model in that such sources lie in a low probability 
corner — a combination of the drop of the UV LF toward 
high luminosity and the drop of the apparent EW dis- 
tribution function toward high EW value. LAE surveys 
with large volume will test the interpretation. 

Therefore, the observed properties of LAEs can be ex- 
plained by simply invoking Lya radiative transfer: the ef- 
fects of the local IGM environment, depending mainly on 
the gas density and line-of-sight velocity and their line- 
of-sight gradients, lead to the distribution of observed 
Lya emission properties at fixed intrinsic Lya luminosity. 
This environmental selection also causes new features in 
the clustering of LAEs that we will study in Paper II. 

8.2. Implications and Discussion 

Our interpretation of the observations of LAEs does 
not invoke any mass dependent dust absor ption, which 
is in contrast to many previous models (e.g.. lDaval et ahl 
2010). Uniformly distributed dust efficiently absorbs 
Lya photons, since the large number of resonant scat- 
terings increase the path length. There is much less 
attenuation when the dust is in gas cl umps and Lya 
photons bounce off the cloud sur faces (jNeufeldl Hflfllt 
lHansen fc Ohl 120061 ). Optical, UV, and Lya obser- 
vations of local star-forming galaxies provide evidence 
that ISM kinematics and geometry play a more signifi- 
cant role than dust in affecting the Lya emission (e.g ., 
IGiavalisco era!ll^lKeel2l)0a 

Our model successfully reproduces the observed UV LF 
of LAEs by incorporating only a mass independent effec- 
tive UV extinction of at most 0.3 mag. We conclude that 
any mass dependent dust effects are not likely to play a 
substantial role to determine the observed properties of 
LAEs, compared to Lya radiative transfer effects. 

Our model also has important implications for the 
duty cycle and the Lya escape fraction of LAEs. The 
theoretically predicted (intrinsic) Lya LF, which essen- 
tially is the halo mass function, is substantially higher 
than the observed one. Two scenarios have been intro- 
duced to address this problem, the duty cycle and the 
Lya escape fraction scenarios (e.g., IStark eiTaD I2007M 
Na gammc et alJl2008f >. In the duty cycle scenario, LAEs 
are short-lived and a fraction of all galaxies are active 
as LAEs at any given time, lowering the amplitude of 
the Lya LF. In the Lya escape fraction scenario, only a 
fixed fraction of Lya photons escape from the source and 
the overprediction problem is solved by shifting the LF 
towards the low luminosity end. To conserve the number 
density of LAEs of a given sample, the masses of host ha- 
los in the duty cycle scenario would be on average lower 
than those in the escape fraction scenario. As a conse- 
quence, the clustering of LAEs would be different in the 
two scenarios, with a stronger clustering in the escape 



fraction scenario. iNaeamine et all (|2008f) find that LAE 
clustering measurements from observations are in favor 
of their duty cycle scenario. 

In our model, Lya photons all escape after a large num- 
ber of scatterings. The Lya escape fraction, in its literal 
meaning, is therefore unity. However, only the central 
part of the extended Lya emission of LAEs can be ob- 
served, which gives rise to an apparent or effective Lya 
escape fraction. Since the observed Lya luminosity has 
a broad distribution at a fixed intrinsic Lya luminosity, 
our model predicts a broad distribution of the effective 
Lya escape fraction rather than a single value. In our 
model, no duty cycle parameter is introduced. Since ha- 
los of the same mass have similar SFR in our model, the 
corresponding intrinsic Lya luminosities are the same, 
i.e., Lya emission does not come from a fraction of ha- 
los. However, an apparent or effective duty cycle arises 
as a result of the selection effect caused by Lya radiative 
transfer (a broad distribution of observed Lya luminos- 
ity at a fixed intrinsic Lya luminosity) and a Lya lumi- 
nosity threshold in observation. This can be seen from 
comparing the UV LF for all galaxies (dropout galax- 
ies and LAEs) and that for LAEs (Fig. [T6|) . which can 
be described as that at a fixed UV luminosity (or halo 
mass) only a fraction of all the galaxies are observed as 
LAEs. This effective duty cycle does not have the phys- 
ical meaning in its original form. Moreover, it is not a 
constant, s i nce it changes with UV luminosity (Fig. [TSJ) . 
iTilvi et all (|2009( ) present an LAE model in which Lya 
luminosity (or SFR) is related to the halo mass accre- 
tion rate, rather than halo mass, and the model naturally 
gives rise to the duty cycle of LAEs. The duty cycle in 
their model, however, has its original meaning, in direct 
contrast with our model. Our model still ties the intrinsic 
Lya luminosity (SFR) to halo mass and let the Lya ra- 
diative transfer do the work of converting it to observed 
Lya luminosity. Because of the large scatter between the 
observed Lya luminosity and the intrinsic one (or halo 
mass) in our model, for a sample of LAEs above a Lya 
luminosity threshold, some of them can reside in halos 
with mass smaller than the threshold mass above which 
halos have the same number density as LAEs. So the ef- 
fective duty cycle in our model has the effect of lowering 
the clustering amplitude of LAEs. 

The main uncertainties in our model are the stellar 
IMF, the SFR, and the intrinsic Lya line profile. The 
first two are general uncertainties for any model. The 
IMF at high-z is neither well constrained observationally 
nor well understood theoretically. The SFR in galaxy for- 
mation model is related to the complex gas physics that 
we do not have a satisfactory understanding. Changing 
IMF or SFR would change the details in the reioniza- 
tion process and therefore change the gas properties (e.g., 
neutral fraction and temperature distribution) at a fixed 
rcdshift. For the reionization history itself, the escape 
fraction of ionizing photons adds a further uncertainty. 
Even though we focus on LAEs at z ~ 5.7, when reion- 
ization is almost complete, different reionization histories 
can still leave different imprints on the gas distribution. 
For example, the IGM temperature in a region is corre- 
lated to the time when this region is reionized and heated 
(e.g., iTrac et aLI 120081 ). A detailed study is needed to 
investigate the effect of inhomogeneous IGM tempera- 
ture distribution on Lya radiative transfer. To be fully 
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self-consistent, for any change in the IMF, SFR, and es- 
cape fraction of ionizing photons, one has to re-run the 
reionization simulation to solve the density, velocity, and 
temperature distributions of neutral gas and then per- 
form the Lya radiative transfer calculation. If the IMF, 
SFR, and escape fraction of ionizing photons change in a 
way to maintain the same reionization history, the effect 
of the IMF and SFR change can be largely character- 
ized by an overall scaling in UV or intrinsic Lya lumi- 
nosity and one does not need to redo the Lya radiative 
transfer calculation. For simplicity and to avoid exten- 
sive computations for reionization and radiative trans- 
fer simulation, we adopt such a scenario in this paper. 
Changing the width of the intrinsic Lya line profile leads 
to changes in the distribution of the apparent to intrin- 
sic Lya luminosity ratio. Although we cast the effect 
as an overall scaling in the apparent Lya luminosity for 
the Lya LF, the intrinsic Lya line profile is important 
in many aspects of the Lya observation (image, spec- 
tra, etc) and its effect deserves a detailed investigation. 
Lya radiative transfer calculation with high resolution 
hydrodynamic simulations for individual LAEs are nec- 
essary to shed light on the intrin sic Lya line profile (e.g., 
lLaursen fc Sommer-Lar scn 2007). Changing the IMF af- 
fects the flux of ionizing photons more than that of the ~ 
1500AUV photons, which leads to a change in the ratio 
of the intrinsic Lya luminosity to UV luminosity, or the 
intrinsic Lya EW. Intrinsic Lya line profile and dust play 
a role in converting the intrinsic EW to apparent EW, 
with the former affecting the Lya luminosity and the lat- 
ter adding extinction to UV luminosity. Therefore, the 
full distribution of observed Lya EW and UV luminosity 
of LAEs can give constraints on the IMF, SFR, dust, and 
intrinsic Lya profile. 

Although we have identified key environment factors 
in shaping the observational properties of LAEs, the de- 
pendence of radiative transfer on environments deserves 
a further study to understand the details of Lya scatter- 
ings in the surrounding regions of LAEs. 

Our Lya radiative transfer calculation relies on the 
gas distribution and properties from the cosmological 
reionization simulation. The radiative transfer of ion- 
izing photons with a ray tracing algorithm is crucial in 
determining the state of gas. We have tested our LAE 
model for a reionization simulation with improved ray 
tracing algorithm (Trac et al. in prep.) in a small box 



(25/i _1 Mpc on a side). We find that the results presented 
in this paper are robust. 

Lya radiative transfer through the surrounding cir- 
cumgalactic and intergalactic media is a physical process 
that likely plays an important role in galaxies at all red- 
shifts. It has to be taken into account for modeling LAEs 
and for interpreting observations. Our model is rather 
simple and can naturally explain an array of observa- 
tions of LAEs, which make it extremely attractive. It is 
interesting to see how well it does in interpreting obser- 
vations of LAEs at lower redshifts (e.g., z ~ 3). We also 
plan to apply it to the era of the late stage of reionization 
to study how to use LAEs to constrain reionization. 
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APPENDIX 

EFFECTS OF THE INITIAL Lyo LINE WIDTH AND GAS TEMPERATURE IN THE exp(-T„) MODEL 

In §[5j we make comparisons of the results from full radiative transfer calculation and those from the exp(— t„) model 
(see Figures IT01 [T3"|) . The exp(— t„) model has been widely adopted in studying LAEs. While the exp(— r„) model we 
present in this paper predicts that the apparent Lya luminosity is at the level of a few percent of the intrinsic one, 
leading to large offset of the apparent Lya luminosity funct i on with respect to the one measured from observation, 
the exp(— Ty) model in other work (e.g.. iMcQuinn et"aLll2007Hlliev et al.ll2008| ) does not show such a large suppression 
in Lya flux. For ex ample, the apparent-to-intrinsic Lya luminosity ratio is at the level of tens of percent for models 
in llliev eT al. (2008j). The appendix aims at resolving the discrepancy by clarifying the assumptions adopted in our 
cxp(— t„) model and others. 

In our exp(— t„), we assume the same initial Lya line width as with the full radiative transfer calculation. The 
initial rms line width is determined by the halo virial temperature (see § [3|) , o"i n j t = 31.9[M fl /(lO lo /i- 1 M )] 1 / 3 kms- 1 . 
We also adopt the gas temperature provided by the reionization simulation, which is solved self-consistently in the 
simulation. As shown in Figure |2T)1 the temperature of the gas drops from ~ 10 5 K near the source center to ~ 10 4 K 
in the IGM. For comparison, the i nitial line width a nd ga s temperature adopted in other work are different from ours. 
For example, the fiducial model in IMcQuinn et al.l (|2007l ) adopt a line width set by the circular velocity at halo virial 
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Fig. 22. — Effects of the initial Lya line width and gas temperature in the exp(— t„) Model. Left panel: Ratio of the apparent (observed) 
to intrinsic Lya luminosity as a function of the initial Lya line width <7i n n. The bottom axis marks the line width in terms of the width <r v ir 
set by halo virial temperature, while the top axis in units of the circular velocity v c at halo virial radius. The thick and thin solid curves 
are the median and quartiles of the luminosity ratio distribution for the case using gas temperature in the reionization simulation. The set 
of dotted curves are the same but with the gas temperature set to 10 4 K in all regions. Right panel: Thermally averaged cross-section of 
Lya scattering as a function of gas temperature. The bottom axis marks the frequency offset with respect to the Lya line center and the 
top axis indicates the fluid velocity that can give rise to such an offset. The cross-section is composed of a thermal (Gaussian) core and a 
Lorentzian wing dominated at small and large frequency offsets, respectively. Changing gas temperature from 10 5 K to 10 4 K has a large 
effect on the Lya scattering cross-section for gas with velocity in the range of ~10-100km s — 1 . See text for details. 



radius, which is about 2.3 times higher than our value. Illiev et al.l (2008) assume the initial rms line width to be 
160km s _1 for all sources, which is much larger than the value we use for most sources. Both lMcQuinn et al.l ()2007l ) 
and Illiev et al.l (|2008|) assume the gas temperature in all regions to be 10 4 K when computing the transmission of Lya 
emission. Given these different assumptions used in our work and others, we perform tests of the effects of initial line 
width and gas temperature on the Lya flux suppression in the exp(— t„) model with a small simulation (box size of 
25/i _1 Mpc on a side). We calculate the distribution of the apparent-to-intrinsic Lya luminosity ratio for all sources 
with halo mass above 5 x 10 9 /i _1 M Q . 

The solid curves in the left panel of Figure [22] show the median and quartiles of the apparent-to-intrinsic Lya 
luminosity ratio as a function of the initial line width for the case using gas temperature in the reionization simulation. 
The results presented in § [5] are based on the (default) value <7i n it /<T V i r = 1. As expected, the luminosity ratio is 
sensitive to the initial line width. At Cinit / 0Vir = 2 the median ratio is about one order of magnitude higher than that 
from our default case. 

The dotted curves in the left panel of Figurel2"2"lcorrespond to the case with gas temperature set to 10 4 K in all regions. 
Changing the gas temperature leads to changes in the thermally averaged cross-section for Lya scattering. As shown 
in the right panel of Figure [22l the scattering cross-section is composed of a thermal (Gaussian) core and a Lorentzian 
wing dominated at small and large frequency offsets, respectively. The core is wider at higher temperature. Since the 
temperature around a source in the reionization simulation can be higher than 10 4 K (Fig. 120 j). setting the temperature 
to 10 4 K would underestimate the scattering optical depth near the source and overestimate the transmitted flux. 
This explains why the luminosity ratio in the case with T = 10 4 K is highe r than that with simulation temperature 
(left panel). For a quantitative comparison, we note that Illiev et a l. (2008) present a case with a- ln it /<7Vir — 2.2 and 
T = 10 4 K (rightmost panel of their Fig. 25) and the median luminosity ratio is ^0.33. Under the same assumptions, 
our exp(— t„) model in Figure |2"21 gives a median ratio of ~0.25, in broad agreement with their result. The residual 
difference may be caused by the difference in the simulations (e.g., hydrogen neutral fraction). 

To conclude, the outcome of the exp(— r„) model depends on the initial Lya line width and gas temperature 
distribution and our exp(— t„) model results are consistent with others if the same assumptions are adopted. 
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